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ABSTRACT 

The  initial  phase  in  a  systematic  analysis  of  weapon- 
ammunition  interaction  has  been  accomplished  under  the  guidance 
of  the  Research  Directorate,  Weapons  Laboratory  at  Rock  Island. 
Acceptance-test  data  for  five  manufacturers'  production  of  5.56mm 
ammunition  were  analyzed  through  time- series  modeling,  an 
empirical  cumulative  distribution  function  was  formulated,  and 
a  bivariate  histogram  of  chamber  pressure  and  port  pressure  was 
developed  for  use  in  the  selection  of  weapon- test  ammunition. 

Statistical  experimental  design  procedures  based  on  facto¬ 
rial  or  fractional  factorial  approaches  are  included  for  use  in 
tests  to  identify  the  controlling  parameters  in  weapon-ammunition 
interactions  and  to  determine  whether  these  parameters  can  be 
identified  from  ammunition  acceptance- test  data.  Some  preliminary 
correlation  analyses  are  included  for  comp-  rison  of  pressure 
measurements  by  means  of  crusher  gages  and  piezoelectric  gages. 


FOREWORD 


This  report  was  prepared  under  Contract  DAAF37060073  by 
Professor  S.  M.  Wu,  Department  of  Mechanical  Engineering  and 
Department  of  Statistics,  University  of  Wisconsin,  Madison,  under 
guidance  of  the  Research  Directorate,  Weapons  Laboratory  at 
Rock  Island,  with  W.  L.  Dahl  as  Project  Scientist. 

The  U.  S.  Army  Small  Arms  Systems  Agency  supported  the  work 
as  part  of  the  Army  Small  Arms  Program  task  entitled  ’'Define 
Weapon  Factors  in  the  Broad  Spectrum  of  Ammunition.” 

The  testing  of  developmental  weapons  has  often  been  limited 
to  the  firing  of  available  quantities  of  ammunition  for  which 
little  information  may  be  available  regarding  either  the  specific 
characteristics  of  the  test  ammunition  or  its  relationship  within 
the  spectrum  of  production  ammunition. 

Definition  of  the  ammunition  characteristics  that  are  im¬ 
portant  in  the  operation  of  differing  gun  mechanisms  and  identi¬ 
fication  of  the  ranges  of  the  variables  will  permit  more  meaning¬ 
ful  tests  of  the  weapons.  These  tests  will,  in  turn,  enable  the 
establishment  of  more  realistic  boundary  conditions  for  analysis 
by  parameter  variation  with  computerized  mathematical  models. 

Such  improved  analytical  techniques  will  provide  better  pre¬ 
dictions  of  weapon-ar.munition  interaction  and  better  estimates 
of  system  reliability  during  the  earliest  stages  cf  the  develop¬ 
ment  process. 


Results  of  the  first  phase  in  an  analysis  of  5.56mm  weapon- 
ammunition  interaction  are  described  in  this  report..  This  por¬ 
tion  of  the  work  has  been  focused  on  the  development  of  methods 
for  the  identification  of  groups  of  sequentially  manufactured 
ammunition-production- lots  that  apparently  fall  within  single 
distributions  of  statistically  normal  production.  The  unique 
control  limits,  associated  with  these  distributions ,  are  being 
explored  as  a  basis  for  the  selection  of  gun- test  ammunition. 

In  addition,  some  statistical  experimental  des.gn  techniques 
are  outlined  for  use  in  the  next  phase  of  the  work.  The  forth¬ 
coming  effort  will  determine  (through  coordinated  tests  at  Frank- 
ford  Arsenal  and  Rock  Island)  whether  the  ammunition  acceptance- 
test  crusher-pressure  data  on  which  this  report  is  based  are 
sufficiently  difinitive  for  determination  of  the  "gun-powering" 
characteristics  of  the  ammunition.  Toward  the  latter  objective, 
some  preliminary  correlation  analyses  of  pressure  data  acquired 
by  means  of  piezoelectric  and  crusher  gages  are  included  in  this 
report. 

One  of  the  appendices  contains  a  description  of  a  preliminary 
mathematical  model,  "Empirical=Mechanistic  Model  for  Interior 
Ballistics  of  Guns,"  formulated  through  cooperation  between  the 
University  of  Wisconsin  and  the  Badger  Army  Ammunition  Plant. 
Although  that  model  was  not  developed  under  the  Research  Directorate 
contract,  it  was  contributed  by  the  University  as  a  possible 
means  by  which  the  technology  may  be  advanced. 
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1.  INTRODUCTION 

Ammunition  and  weapon  need  to  be  compatible.  Ammunition  tests 
are  performed  to  ascertain  that  the  ballistic  characteristics  of  ammunition 
are  consistent  with  the  gun  design  requirements.  The  tests  are  useful 
only  if  the  responses  provide  meaning  al  information  toward  this  end.  The 
present  methods  of  testing,  namely  copper  crusher  and  piezo  methods, 
need  to  be  compared  to  determine  the  extent  of  useful  information  obtained . 

Copper  crusher  forms  the  present  method  of  acceptance  testing. 

The  control  limits  used  for  acceptance  are  based  upon  specifications  which 
may  or  may  not  represent  the  real  capabilities  of  the  present  production 
process.  Additionally,  the  control  limits  may  be  only  indirectly  related  to 
the  functional  requirements  of  automatic  weapons.  It  is  therefore  possible, 
that  unnecessarily  severe  requirements  for  versatility  might  have  been  im¬ 
posed  on  the  weapon.  Analysis  of  acceptance  testing  data  is  required  to 
determine  more  realistic  control  limits. 

Just  as  ammunition  tests  are  performed  to  evaluate  ammunition 
characteristics,  similarly  weapon  tests  are  needed  to  evaluate  weapon  per¬ 
formance,  which  must  lie  within  specified  limits.  Since  the  ammunition  lots 
are  not  identical  in  their  ballistic  properties,  a  criterion  has  to  be  established 
to  select  ammunition  lots  for  weapon  testing. 

The  first  step  toward  the  establishment  of  such  a  criterion  is  to  de¬ 
termine  the  ammunition  characteristics  that  influence  weapon  performance. 
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The  capability  of  ammunition  manufacturing  process  can  then  be  analyzed 
in  terms  of  these  characteristics  to  evaluate  the  differences  between 
ammunition  lots .  Lots  that  would  give  a  large  variation  in  weapon  perfor¬ 
mance  would  be  selected  for  weapon  testing. 

Related  to  the  analysis  of  ammunition  manufacturing  process  and 
weapon  testing  are  the  questions  of  determining  proper  number  of  tests  to 
be  conducted  at  each  stage  of  data  generation  and  detection  of  changes  in 
ammunition  characteristics  at  different  stages  of  manufacture .  The  latter 
would  help  explain  the  final  ballistic  characteristics  attained  by  the  am¬ 
munition.  It  is  felt  that  sufficient  attention  has  not  yet  been  give  to  these 
questions,  satisfactory  solutions  of  which  are  ilkely  to  lead  to  considerable 
savings  to  the  government. 

It  is  felt  that  a  theoretical  understanding  of  the  internal  ballistic 
mechanism  involved  would  be  very  useful  in  the  interpretation  of  the  results 


obtained . 
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SUMMARY 

The  analyses  in  this  report  are  based  upon  acceptance 
test  data  for  5.56  mm.  Ball  M-193  and  5.56  mm.  Tracer  M-196 
ammunition  from  five  manufacturers  (Lake  City,  Twin  City, 
Remington,  Federal  and  Winchester)  covering  a  period  from 
Jply  1968  to  March  1971.  Additionally,  special  ammunition 
test  data  from  Frankford  and  ammunition  data  from  B.A.A.P. 
have  bee:,  analyzed.  Acceptance  test  data  from  B.A.A.P. 
covers  a  production  period  of  the  past  three  years.  Copper 
crusher  and  piezo  data  for  the  same  propellant  lot  are  also 
available. 

The  purposes  of  these  analyses  are  many  fold.  They 
include  a  comparison  of  copper  crusher  and  piezo  methods  of 
testing,  the  determination  of  more  realistic  control  limits 
for  ammunition  production  and  selection  of  ammunition  lots 
for  weapon  testing.  Chapter  1  brings  out  the  importance  of 
such  analyses. 

Piezo  transducer  gives  more  information  regarding  pres¬ 
sures  inside  the  barrel.  Maximum  pressure  at  the  gage  loca¬ 
tion  as  well  as  ignition  delay  and  slope  can  be  determined. 
These  are  indicators  of  propellant  characteristics.  Crusher 
deformation  can  be  considered  as  a  weighted  integral  of 
piezo  pressure-time  curve  and  could  be  a  good  estimate  of 
impact  energy.  Analysis  in  Chapter  2  points  out  that  copper 
crusher  indicates  similar  pressure-velocity  relationship 
as  piezo  does  (Figs.  7  and  10),  but  it  has  larger  variability. 
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Coating  data  against  copper  crusher  chamber  pressure  plot 
does  not  show  any  pattern  (Fig,  8}  ,  whereas  plot  of  piezo 
peak  chamber  pressure  versus  coating  date  reveals  a  time 
trend  (Fig.  11).  The  time  trend  would  be  lost  if  only  cop¬ 
per  crusher  data  are  examined. 

Chapter  3  contains  a  detailed  analysis  of  standard 
test  data.  Chamber  pressure  data  from  standard  tests  are 
found  to  exhibit  marked  trends.  Three  methods  (semi  aver¬ 
ages,  cumulative  sum  and  moving  averages)  have  been  employed 
to  estimate  the  underlying  process  behavior.  In  particular, 
the  point  of  shift  after  which  the  trend  is  less  predominant 
is  determined.  The  method  of  cumulative  sum  is  found  to 
give  the  best  visual  indication  of  the  point  of  shift. 

The  chamber  pressure  data  are  found  to  be  nonstationary 
(Lake  City  data)  even  after  the  point  of  shift.  The  data 
are  also  found  to  be  autocorrelated .  Therefore,  time  series 
models  have  been  obtained  to  explain  the  nature  of  correla¬ 
tions.  Analysis  of  chamber  pressure  data,  from  different 
manufacturers  shows  the  mean  chamber  pressures  to  be  quite 

close  to  each  other.  The  chamber  pressure  variance  is  found 

4  2 

to  vary  considerably  (118  x  10  psi  for  Remington  to  629  x 
10^  psi2  for  Winchester)  from  manufacturer  to  manufacturer. 
The  mean  chamber  pressure  of  Ball  ammunition  (46600  psi)  is 
lower  than  that  of  Tracer  ammunition  (49200  psi) . 

Analysis  of  chamber  pressure  variance  shows  a  larger 

5  2 

value  of  within  lot  variance  (29  x  10  psi  )  for  ammunition 
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from  Twin  City  as  compared  to  (17.6  x  10  psi  )  ammunition 
from  Lake  City.  In  both  cases  the  trend  is  towards  a 
reduction  in  variance,  indicating  continued  improvement  in 
production  and  testing  processes. 

In  Chapter  4,  a  method  based  upon  empirical  cumulative 
|  *  distribution  function  has  been  developed  to  obtain  control 

limits  for  ammunition  production  processes.  The  control 
limits  are  based  upon  99  percentile  point  of  the  empirical 
cumulative  distribution  function.  To  narrow  the  control 
limits,  data  after  'cutoff  date'  alone  have  been  considered. 
The  conditions  under  which  cutoff  date  can  be  taken  as  the 
date  corresponding  to  the  point  of  shift,  are  given  in 
Section  4.2.  Cutoff  date  and  control  limits  have  been  cal¬ 
culated  for  Ball  and  Tracer  chamber  pressure  data  from  the 
five  manufacturers.  On  the  average  the  control  limits  have 
been  reduced  by  2000  psi  from  the  existing  ones. 

In  Chapter  5,  a  criterion  has  been  developed  for  select¬ 
ing  ammunition  lots  for  weapon  testing.  It  is  based  upon 
the  bivariate  histogram  of  chamber  pressure  and  port  pressure 
from  standard  tests.  The  inadequacy  of  basing  the  selection 
on  chamber  pressure  or  port  pressure  alone  has  been  discussed. 
The  selected  lots  have  been  classified  into  High,  Medium  and 
Low.  Medium  lots  are  based  upon  the  mode  and  the  High  and 
the  Low  lots  are  based  upon  the  approximate  90%  confidence 
limits  of  the  bivariate  histogram.  An  experimental  procedure, 
using  fractional  factorial  or  factorial  designs,  has  been 
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suggested  to  determine  ammunition  parameters  that  control 
weapon  performance. 

Chapter  6  contains  a  suggested  procedure  to  determine 
the  number  of  tests  necessary  to  estimate  ammunition  para- 
JHete^s,  The  procedure  is  based  upon  the  desired  precision 
of  estimates  and  the  experimental  error  involved.  The  se¬ 
quential  method  for  estimation  of  experimental  error  shows 
that  twenty  readings  do  not  give  a  proper  estimate  of  ex¬ 
perimental  error  in  copper  crusher  testing  and  ten  readings 
are  not  sufficient  to  obtain  a  good  estimate  of  experimental 
error  in  piezo  testing.  The  current  practice  in  standard 
testing  is  to  obtain  an  estimate  of  standard  deviation  based 
upon  20  tests.  The  implication  of  this  analysis  is  that  the 
estimate  so  obtained  is  likely  to  be  modified  considerably 
if  it  is  based  on  sufficiently  large  number  of  tests. 

Chapter  7  deals  with  the  analysis  of  data  at  different 
stages  of  manufacture.  Analysis  of  propellant  lot  charac¬ 
teristics  from  B.A.A.P.  indicates  the  lot  characteristics 
(charge  weight  and  chamber  pressure)  to  be  serially  correlated. 
Comparison  of  acceptance  test  results,  for  the  same  propellant 
lots  from  B.A.A.P.,  Lake  City,  and  Twin  City  show  the  test 
results  to  be  almost  identical,  indicating  no  change  in 
propellant  properties  during  the  time  period  between  the 
tests. 

During  the  course  of  the  project,  need  was  felt  for  a 
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[  '  model  for  the  interior  ballistic  system  of  guns.  Appendix 

! 

|  I  contains  an  initial  attempt  toward  building  an  empirical- 

i 

|  mechanistic  model,  using  a  Lagrangean  formulation  of  the 

f 

| 

i  hydrodynamic  system.  Piezo  pres sur e-time  curve  and  velocity 

I  . 

|  have  been  used  as  responses.  The  model  is  shown  capable  of 

i  iterative  improvements.  Several  possible  uses  for  the  model 

|  have  been  outlined. 

The  data  used  in  different  analysis  in  this  report  are 
i  given  in  Appendices  II,  III  and  IV.  Computer  programs  neces- 

|  1  sary  for  the  analyses  are  included  in  Appendix  V. 
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2.  COMPARISON  OF  COPPER  CRUSHER  AND 
.  PIEZO  METHODS  OF.  TESTING 


^  '  *«  'j  vVv  *  ' 

OooDer  crusher  and  Diezo  methods 


ammunition  tests .  The  crusher  gage  has  been  in  vogue  for  almost  a 
century  and  is  in  use  even  today  as  the  sole  standard  method  of  pressure 
measurement.  The  piezo  gage,  even  though  known  for  a  considerable  period 
of  time,  is  still  not  adopted  as  a  standard  method  of  pressure  measurement. 
As  a  result,  relatively  small  amounts  of  piezo  data  is  available  compared 
to  the  large  amount  of  copper  crusher  data  accumulated  over  past  years. 

In  this  section,  the  two  methods  are  compared  regarding  their  relative  use¬ 
fulness. 
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2.1.  Purpose  of  Testing 

Purpose  of  testing  is  two  fold:  as  a  means  of  acceptance  testing 
and  as  a  tool  for  process  control.  Comparison  is  piimarily  based  upon  the 
former  function  of  the  testing  method.  Where  an  ammunition  lot  is  tested 
for  acceptance,  it  is  necessary  to  determine  whether  the  powder  can  impart 
desired  velocity  to  the  bullet  and  whether  the  gun  can  withstand  the  pres¬ 


sures  generated. 


2.2  Description  of  Piezo  Responses 

Piezo  response  is  acontinuous  pressure-time  curve  for  the  particu¬ 
lar  section  along  the  barrel  (usually  mid  chamber  position  for  chamber 

>  •  i  -i  i  -  *  *  ’  i  ***'* '  •  * 

pressure  measurement)  where  it  is  located.  One  typical  curve  is  shown 

c  ;  '  •  •  *  * 

in  Fig.  1.  There  is  an  initial  portion  of  'ignition  delay'  during  which  the 
burning  rate  is  small.  Next  there  is  a  rapid  rate  of  rise  of  pressure  due 
to  increased  burning  rate.  The  expansion  of  gases  behind  the  bullet  has 
a  tendency  to  reduce  the  pressure.  Eventually  the  powder  burns  off  and 
as  a  result  of  these  interacting  causes,  the  pressure  reduces.  The  pres¬ 
sure  curve,  therefore,  exhibits  an  unimodal  maximum. 

Piezo  transducer  has  a  time  constant  of  the  order  of  10"9  seconds 
and  is  quite  widely  used  to  obtain  the  pres  sure -time  history. 

2 . 3  Copper  Crusher  Response  and  its  Relationship  with  Piezo  Response 

In  standard  acceptance  testing  for  measurement  of  chamber  pressure, 
the  pressure  inside  the  gun  chamber  is  transmitted  to  the  copper  cylinder 
through  a  rigid  steel  piston.  The  pressure  acting  on  the  copper  cylinder 
is  the  same  as  in  Fig.  1.  If  the  elastic  and  plastic  behavior  of  the  copper 
cylinder  is  known,  the  crusher  deformation  can  be  expressed  as  a  weighted 
integral  of  the  piezo  pressure-time  curve. 
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2 . 4  Method s  of  Choosing  the  Acceptance  Test  Criteria 


Gun  body  may  fail  if  the  burning  conditions  ,are  severe .  One  way 
to  estimate  the  limiting  condition  is  by  means  of, the  maximum  stress  pro¬ 
duced  in  the  gun  walls  which  is  a  function  of  the  maximum  pressure  devel¬ 
oped.  The  peak  piezo  chamber  pressure  is  perhaps  a  suitable  indicator  ,  i 

1  i 

for  this  type  of  failure.  Another  possibility  is  impact  failure,  which  appears 
more  realistic  for  the  present  situation.  Deformation  of  copper  cylinder  is 
one  of  the  ways  of  measuring  impact  energy.  However,  the  best  indicator 
of  impact  energy  can  only  be  determined  by  a  detailed  solid  mechanical 

I 

analysis  of  the  gun  body.  . 
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2 . 5  Correlation  Analysis 


Correlation  analysis  was  conducted  to  determine  the  type  of  information 
that  can  be  obtained  from  the  data.  Three  major  sets  of  data  were  used,  1 

‘  i 

i 

rarely,  standard  acceptance  test  data,  special  ammunition  test  data  supplied 

i  i 

by  Badger  Army  Ammunition  Plant  (B.A.A,  P.) .  These  arc  listed  in  Appendix  II. 


M 
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2.5 . 1  Standard  Acceptance  Test  Data 


Pertinent  results  in  this  category  are  summarized  below: 

(1)  Chamber  pressure  has  no  correlation  with  velocity  from  velocity 


i<  i 


barrel.  (Fig.  2) 

! 

(2)  Chamber  pressure  and  port  pressure  have  a  slight  negative  correlation 


(Fig.  3). 


Chamber  Pressure 
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2.5.2  Special  Ammunition  Test  Data 

Twenty  pieces  of  special  ammunition  test  data  were  obtained 
from  Frankford  Arsenal.  Limited  information  available  is  as  follows: 

.  ,  .(1)  There  is  a  positive  correlation  between  peak  chamber  pressure 

T  ’  '  v 

$  ■* 

and  velocity.  (Fig.  4) 

(2)  Peak  chamber  pressure  is  uncorrelated  with  peak  port  pressure. 

(Fig.  5) 

(3)  Pressure-time  integral  is  uncorrelated  with  velocity.  (Figs.  6(a), 

6(b)  ) 

2.5.3  B.A.A.P.  Data 

Copper  crusher  and  piezo  transducer  data  for  composites  and  hand- 
blends  were  supplied  by  B.A.A.P.  Results  are  summarized  below: 

(a)  Copper  Crusher  Method 

(1)  Chamber  pressure  and  velocity  from  pressure  barrel  are  correlated. 
(Fig.  7) 

(2)  Plot  of  coating  date  against  chamber  pressure  reveals  no  trend . 
(Fig.  8) 

(3)  Plot  of  coating  date  against  velocity  from  pressure  barrel  has  a 
slight  trend .  (Fig .  9) 


Pressure-tiroe  Integral 


.  Coating  Date 
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(b)  Piezo  Method 

(1)  Peak  chamber  pressure  and  velocity  are  correlated.  (Fig.  10) 

(2)  Plot  of  coating  date  against  peak  chamber  pressure  reveals 
a  time  trend.  (Fig.  11) 

0  (3)  Slope  ofspressure-tlme  curve  is  negatively  correlated  with 

ignition  delay  (Fig.  12)  and  is  positively  correlated  with  veloc¬ 
ity  (Fig.  13)  and  peak  pressure.  (Fig.  14) 
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Similar  correlative  structure  is  observed  in  the  case  of  Hand 


Blends. 

1  <  ,  . 
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2. 6  Evaluation  of  Copper  Crusher 

Copper  crusher  has  been  found  to  be  a  reliable  method  of  competitive 
measurement.  It  is  possible  that  the  crusher  deformation  is  a  good  estimate 
of.  impact  energy.  From  Figures  7  and  10  it  can  be  seen  that  the  crusher 
.Indicates  the  same  trend  as  piezo  does,  but  indicates  a  larger  variability. 
Coating  date  against  copper  crusher  chamber  pressure  plot  does  not  show 
any  pattern,  whereas  plot  of  coating  date  versus  piezo  peak  chamber  pressure 
reveals  a  time  trend .  This  time  trend  could  be  a  significant  factor  in  pro¬ 


duction  process  control  and  would  be  lost  if  only  copper  crusher  data  are 


examined . 
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3 .  ANALYSIS  OF  ACCEPTANCE  TEST 'PAT A'  'o:  -j 


‘Acceptance  test  data' for  5 -5  6mm,  Bill  M-193  ammunition  from 
five  manufacturers  (Lake  City,  Twin  City,  Remington,  Federal  and' 
Winchester)  as  well  as  for  5.56mm.  Tracer  M-196  ammunition  from  three 
manufacturers  (Lake  City,  Twin  City  and  Winchester)  were  made  available 
by  Rock  Island  Arsenal  for  investigation.  The  data  were  subsequently  up¬ 
dated  to  cover  a  production  period  from  July  1968  to  March  1971.  The 
acceptance  test  data  contain  the  following  information: 

(1)  Ammunition  lot  number  and  date  testing 

(2)  Propellant  lot  number 

(3)  Average  charge  weight  used  for  the  ammunition  lot 

(4)  Chamber  pressure  -  which  is  the  average  of  20  chamber  pressure  read 
ings  per  ammunition  lot 

(5)  Maximum  value  of  these  20  readings 

(6)  Port  pressure  -  which  is  the  average  of  20  port  pressure  readings  per 
ammunition  lot 

(7)  Muzzle  Velocity  -  which  is  the  average  of  20  velocity  readings  per 
ammunition  lot 

(8)  Standard  deviation  for  chamber  pressure,  port  pressure  and  velocity. 
This  is  computed  from-the  corresponding  20  individual  readings. 

(9)  Correction  for  chamber  pressure  and  velocity 


(10)  Accuracy  -  which  is  a  measure  of  how  accurate  the  flight  of  the 
bullet  is 

(11)  The  data  also  contain  information  regarding  the  reference  lot  used, 
number  of  bullets  fired  using  the  test  pressure  barrel  and  velocity 
barrel,  reference  velocity  and  reference  chamber  pressure  as  well  as 
the  press  're  barrel  number  and  velocity  barrel  number 

3,1  Preliminary  Study  of  the  Data 

Graphical  method  was  used  to  obtain  a  visual  picture  -of  the  data . 
Figures  IS  to  22  show  the  plots  of  chamber  pressure,  port  pressure  and  muz¬ 
zle  velocity  against  the  date  of  testing.  The  plots  include  Ball  ammunition 
data  from  five  manufacturers  and  Tracer  ammunition  data  from  three  main 
Chamber  factories.  Pressure,  port  pressure,  and  muzzle  velocity  were  also 
plotted  against  the  corresponding  ammunition  lot  numbers  and  the  resulting 
plots  are  given  in  Figures  22  to  27. 

A  visual  examination  of  these  graphs  reveals  the  following: 

(1)  Testing  is  not  done  «t  regular  time  intervals.  It  appears,  therefore, 
that  plotting  based  upon  lot  numbers  provides  a  better  representation 
of  the  time  sequence  of  the  production  process. 

(2)  The  data  show  production  trends,  in  particular,  .he  chamber  pressure 
data  from  Lake  City  and  from  Twin  City  show  a  marked  shift. 
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FIGURE  24 
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FIGURE  27 


(3)  In*  addition?to  the  differences^in'ths  trend.there  appear  to; be* dif¬ 

ferences  in  the  range  {range  =  maximum -minimum)  of  ithe  data 
.from  the-five  manufacturers ;  *  •  •>.-->  j>  ■: ■ 

(4)  Mean  chamber  pressure  for  Ball  ammunition  is  lower  than  that  for 
Tracer  ammunition  while  muzzle  velocity  for  Ball  ammunition  is 
higher  than  that  of  the  Tracer  ammunition. 

A  quantitative  analysis  of  the  standard  testing  data  is  now  con¬ 
ducted  . 

3 . 2  Determination  of  the  Point  of  Shift 

It  has  been  noted  that  the  chamber  pressure  data  reveal  the  existence 
of  a  shift,  after  which  the  data  assume  a  relatively  stationary  pattern.  The 
large  variability  in  the  observed  data  makes  it  difficult  to  determine  the 
point  where  the  underlying  process  shifts.  Methods  are  therefore  needed  to 
produce  a  visual  picture  reflecting  the  underlying  process  and  minimize  the 
random  fluctuations  about  it ,  Any  such  method  would  depend  upon  certain 
assumptions  regarding  the  process  and  would  be  only  good  within  the  assump¬ 
tions  made.  Three  methods  of  estimating  the  underlying  process  would  now  be 
considered.  These  methods  would  also  serve  to  determine  a  unique  answer  for 
the  point  of  shift  rather  than  differing  answers  that  would  result  by  eyeball 


estimation . 
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of  the  group  average,  the  resulting  plot  presents  a  better  visual  picture 

v-'/'j  ;r  ; 

of  the  underlying  process.  Ifrom  Figures  28  and  29  it  can  be  seen  that 
large  value  of  results  in  greater  smoothenlng  and  a  better  visual,  pic¬ 
ture.  However,  excessive  smoothenlng  might  lead  to  a  loss  of  information. 
An  optimum  value  of  K^  Is,  therefore,  necessary.  Point  of  shift  is  the  lot 
number  after  which  the  group  averages  stabilize. 

3.2.2  The  Method  of  Cumulative  Sum 

This  is  a  plot  of  successive  partial  sums 


m 


i  =  1 


(Xj  -  KJ  ;  m  =  1,  2— -n 


plotted  against 'm\  Here 

X,  =  chamber  pressure  of  the  ith  lot 
n  =  total  number  of  lots  considered 
and  Kg  =  reference  constant 

The  shape  of  the  plot  changes  with  the  different  values  of  constant 
Kg.  Two  plots  of  the  cumulative  sum  (cusum)  are  shown  in  Figures  30  and  31. 
The  mean  level  over  any  portion  of  the  cusum  plot  is  given  by 


. . Change  in  the  cumulative  sum 

Mean  Level  =  K.  +  - d — - : — r; - 

2  Change  in  M 


In  particular,  an  estimate  of  the  mean  chamber  pressure  for  a  particular  lot 


uw 


GROUP  NO. 


is  given  by  the  slope  of  cusum  plot  at  that  lot  plus  Kg.  A  change  in  slope/ 
therefore,  signifies  a  change  in  the  process  level  and  forms  a  criterion  for 
the  determination  of  the  point  of  shift.  An  optimum  value  of  Kg  is  one  that 
gives  the  best  Indication  of  change  lr.  process  level. 

3.2.3  Method  of  Moving  Average  s 

Moving  average  is  another  technique  used  to  produce  an  approximation 
to  the  underlying  process .  The  plot  is  obtained  as  follows:  An  average  of 
first  Kg  observations  is  computed  (Kg  =  period)  and  is  plotted  at  the  mid- 
period  position.  Now  the  first  observation  is  deleted  and  (Kg  +  1)  th 
observation  is  added  to  get  a  new  average.  This  is  plotted  at  the  mid  point  of 
this  period  and  the  process  is  repeated  till  all  available  observations  are  - 
exhausted . 

Each  point  on  the  curve  is  an  estimate  of  chamber  pressure  for  that 
lot.  Since  plotting  is  done  at  the  mid-point  of  a  period,  Kg  should  be  an 
odd  number  so  that  the  plotted  point  would  correspond  to  an  actual  lot.  Two 
such  plots,  with  Kg  equal  to  5  and  15  are  shown  in  Figures  32  and  33.  It 
can  be  observed  that  the  longer  the  period,  the  more  reduction  in  fluctuation. 
As  the  length  of  the  period  Increases,  there  is  a  tendency  for  the  moving 
average  to  'iron  out'  the  underlying  process.  It  is  prudent,  therefore,  to 
use  as  short  a  period  as  possible  consistent  with  a  reasonably  smooth  moving 


average . 
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f 7 ' "All*  the' thte'e  riietebds  iWnFe^fr  pwTJOife^of  ei&hatihg;tee  under¬ 
lying  process  by  reducing  the  random  fluctuations  about  it.  If  the  sue- 
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cessive  values  of 'chamber  pressure  can  be  assumed 'to  be  uncorrelated, 


then  optimum  values  of 'constants  ;  ’  K^V  Had  caii'be  determined .  For 

example ,  in  the  method  of  moving  averages,  the  underlying  process  can  be 

*  J  ,  • 

estimated' for  any  particular* value  of  K.  .  The  deviation  of  the  observations 

I  '  i  3 

•  I 

from  this  Urtderr  /ing‘ process  can 'be  coiiiputed.  The  value  of  K„  should  be 

•3 

so  chosen  that. the  variance  of  these  deviations  is  equal  to  the  average 

J  * 

Internal  (within  lot)  variance  of  the  mean  chamber  pressure.  This  has  not 

been  actually  carried  out  since  the  chamber  pressure  readings  are  found  to 

!  ’  ? 

i  constitute  a  time  series  of  correlated  observations. 

1  1 

The  purpose  of  determlng  the  point  of  shift  (the  date  on  which  the 

4 

shift  occured  is  termed  the  'cutoff  date')  is  to  be  able  to  consider  only 

i 

■  the  stationary  piart  of  the  data .  Control  limits  based  upon  data  after  the 

!  ;  . 

cutoff  date  are  expected’  to  be  harrower  than  the  present  control  limits .  If 

•  i 

the  data  are  plotted  a  rough  estimate  of  the  point  of  shift  can  be  obtained 
visually.  Referring  to  Figures  28  through  33,  it  is  observed  that  the  point 

,  i 

of  shift  is  best  determined  by  the  cumulative  sum  plot,  A  good  value  for  the 

i 

reference  constant  K„  seems  to  be  the  mean  of  the  entire  series  of  data 
because  this  value  of  reference  constant  gives  a  large  change  in  slope  at 
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the  point  of  shift.  This  procedure  wiU'be'used  ^  4  to  obtain  the 

cutoff, dat^for, the  chamber  pressure  date  for  Ball  and  Tracer  ammunition. 

)  t  4  •  -  s  *  -  \  t 

3.3  Stationarity  of  Lake  City  Chamber  Pressure  Data  After  Cutoff 

i "  v  . 

is  tp.be  determined  if  the  chamber  pressure  series  is  stationary 
after  cutoff.  The  chamber  pressure. data  consist  of  a  time  series  of  obser¬ 
vations  .which  are  means  of  20  chamber  pressure  readings.  Strict  station¬ 
arity  demands  that  all  these  readings  come  from  the.  same  distribution,  in 

‘  v 

particular,  they  have,  the  same  mean  and  the  same  variance. 

...  Lake  City  chamber  pressure  variance  can  be  considered  to  be  approx¬ 
imately  constant  (It  is  not  exactly  constant,  as  will  be  seen  later)  and  has 

2 

a  mean  value  of  1760806  (psi)  .  The  variance  of  mean  chamber  pressure  is, 

2 

therefore,  equal  to  88040  (psi)  (»  1760806/20) .  This  is  a  measure  of  the 
average  internal  (within  a  lot)  variance  of  the  mean  chamber  pressure.  A 
measure  of  variance  between  the  mean  chamber  pressures  of  different  lots  can 
be  obtained  by  calculating  the  variance  of  the  chamber  pressure  series  after 
cutoff.  .This  is  found  to  be  462000  (psi)  .  If  all  the  chamber  pressure  read¬ 
ings  had  the  same  mean,  then  the  variance  between  lots  would  be  approxi¬ 
mately  the  same  as  variance  within  the  lot.  In  reality  this  ratio  (452000/ 
88040)  is  of  the  order  of  50,  hence  chamber  pressure  data  after  cutoff  is  still 
nonstatfenary. 
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3.4  NonNormality  of  Lake  City  Chamber  Pressure  Data 


It  is  of  Interest  to  see  if  the  chamber  pressure  readings  constitute 

observations  from  a  normal  distribution.  Figure  34  shows  the  histogram  for 

Lake  City  chamber  pressure  data.  The  best  possible  normal  distribution  to 

fit  the  data  is  found  to  be  N(462,  12.7  ).  Table  1  shows  the  Goodness  of 

Fit  test  to  determine  if  the  data  fit  this  normal  distribution.  Since  the  cal- 

2 

culated  chi-square  value  of  23.76  is  greater  than  x  ^  3  o  05  (=16.9), 
normality  assumption  Is  not  justified  at  95%  level  of  confidence.  Figures 
35  to  38  show  histograms  of  chamber  pressure  data  (after  cutoff)  from  Twin 
City,  Winchester,  Ramington,  and  Federal.  The  histograms  are  clearly 
seen  to  be  nonnormal.  Because  of  the  presence  of  trend  even  after  cutoff, 
the  underlying  normal  distribution  of  each  observation  is  distorted .  The 
result,  is,  therefore,  not  unexpected . 

3.5  Normality  of  Lake  City  Port  Pressure  Data 

Histogram  of  port  pressure  readings  after  cutoff  is 
shown  in  Fig.  39,  Table  2  shows  the  Goodness  of  Fit  test  to 
determine  if  the  data  fit  the  estimated  normal  distribution. 

The  calculated  chi-square  value  is  113.4  and  is  greater  than 
the  critical  value  18.3  (x^3-3  qq5  *  18.3).  Normality  assump¬ 
tion  is  not  justified  at  95%  level  of  confidence. 
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Figure  34 
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Goodness  of  Fit  Test  for  Lake  City  Chamber  Pressure  Data 
(After  Some  Regrouping) 


Interval 
)e  1 00  psi 


Expected  no.  ‘  ““ 
of  occurances 

Observed  no.  n8i  .2 

of  ocgurances  Fitted  Normals 

£i _ N(462,12.7)2 _ n6^ 


372.0 

to 

430.95 

4 

6.6 

1.00 

430.95 

to 

437.5 

9 

17.9 

4.35 

437.5 

to 

444.05 

48 

51.0 

0.18 

444.05 

to 

450.6 

97 

99.0 

0.04 

450.6 

to 

457.15 

204 

162.0 

10.90 

457.15 

to 

463.7 

186 

188.5 

0.03 

463.7 

to 

470.25 

177 

177.0 

0.00 

470.25 

to 

476.8 

105 

127.0 

3.80 

476.8 

to 

433.35 

69 

70.0 

0.16 

483.35 

to 

489.9 

28 

30.0 

0.13 

489.9 

to 

496.45 

9 

10.3 

0.17 

496.45 

to 

503.0 

6 

3.0 

3.00 

Total= 

=942 

Total=942  Total=> 

23.76 
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Since  the  chamber  pressure  and  port  pressure  data  are  not  normally 
distributed,  an  empirical  approach  has  been  used  for  the  calculation  of 
control  limits  as  well  as  for  ammunition  selection  for  weapon  tests. 

3.6  Comparison  of  Chamber  Pressure  Data  from  Different  Manufacturers 

Table  3  shows  the  mean  and  variance  of  chamber  pressure  data 
from  different  manufacturers .  Data  after  the  cutoff  date  alone  have  been 
used  for  this  comparison.  The  variance  (between  lots)  of  chamber  pressure 
is  seen  to  vary  considerably  from  manufacturer  to  manufacturer.  Means 
are  observed  to  be  quite  close  to  each  other.  The  mean  chamber  pressure 
for  Ball  ammunition  is  lower  than  that  for  Tracer  ammunition. 
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TABLE  2 

Goodness  of  Fit  Test  for  Lake  Pity  Port  Pressure  Data 
(After  Some  Regrouping) 


Interval 
x  10  psi 

Observed  No. 
of  Occurances 
fi 

Expected  No. 
of  Occurances 
n01 

Fitted  Norma  L 
N(1533,  36.6^ 

(Fi-nfr)2 

nOi 

1300  to  1390 

6 

1 

25 

1390  to  1426 

5 

2 

4.5 

1426  to  1444 

8 

5.6 

1.0 

1444  to  1462 

13 

17 

1.0 

1462  to  1480 

35 

45 

2.2 

1480  to  1498 

41 

91.4 

27.8 

1498  to  1516 

146 

146 

0 

1516  to  1534 

217 

164 

17.1 

1534  to  1552 

236 

185 

14.0 

1552  to  1570 

135 

139 

0.1 

1570  to  1588 

44 

84 

19.0 

1588  to  1606 

41 

41 

0 

1606  to  1660 

15 

21 

1.7 

Total  =  942 

Total  =  942 

Total  =113.4 

TABLE  3 


Comparison  of  Mean  and  Variance  of 
Chamber,  Pressure  from  Five  Manufacturers 
(Data  after  Cutoff  Date) 


Manufacturers 

Chamber  Pressure  (Ball) 

Chamber  Pressure  (Tracer) 

Mean 

Variance  2 

Mean 

Variance  2 

x  100  psl 

x  10,  000  psi 

x  100  psi 

x  10,  000  psi 

Lake  City 

462 

462 

488 

321 

Twin  City 

461 

217 

490 

236 

Winchester 

461 

629 

499 

140 

Remington 

474 

118 

Federal 

474 

238 
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3. 7  Time  Series  Analysis  of  Lake  City  Data  after  Cutoff  )  i 

l 

Time  Series  Analysis  was  carried  out  on  Lake  City  Chamber 
pressure  data  after  cutoff.  A  total  of  942, observations  were  : 
considered.  The  model  was  found  to  exhibit  a  nonstationary 
seasonal  behavior.  The  fitted  model  is 

(1  -  4»22B22  r  4>23B23)  (1  -  B)Zt  =  (1  -  61B)at 

> 

1.1. 

where  is  the  chamber  pressure  at  t  lot.  i 

The  estimated  parameter  values  are 

<J>22  =  -0.05613 
J23  =  0.1307 

6X  =  0.6651 

The  parameter  confidence  intervals  are 

-0.1169  <  <f>22  <  0.0046  1 

0.0704  <  4>23  <  0.1910 
0.6161  <  0,  <  0.7142 

A  \  ' 

Even  though  the  confidence  interval  on  4>22  includes  zero, 

a  much  better  fit  is  ob- 


i.OJ 

2 

X  test  indicate  the  model 


it  was  decided  to  retain  it  because 
tained  by  its  inclusion. 

The  correlation  matrix  is 


^22 

1.0 

to 

OJ 

0.0319 

1.0 

0.0309 

0.0312 

The  error  autocorrelations  and 


to  be  adequate. 


e 


.  i 


i  i  . 

Chamber  pressure  variance,  within  each  lot  for  Lake  City  and  for 

.  .  .  ■■  • 

Twin  City  is  plotted  in  Figures  42  and  43^ respectively.  The  average  value 

I  I  ”*  « 

of  within  lot  variance  for  Lake  City  is  1760806-v(psi)  and  for  Twin  City  is 

!  £  1  -  \v 

2  905  842  (psi)  .  Lake  City,  therefore,  has-a-  smaller:within  lot  chamber 

i  1  ..  A 

pressure  variance  as  compared  to  Twin.City.  If  it  were  possible  to  assume 

I  !  ’ 

’  I  ~  v  " 

that  identical  ball  powder  is  supplied  to  both  the  plants,  then  this  difference 

, 

»  .  -  ' 

in  variance  should'  be  attributed  to  variation  in  components  and/or  standard 
testing  procedures  l  ! 

i  . 

The  method  of  cumulative  sum  is  now  employed  to  determine  the  under- 

.  '  1 

lying  process.  The  results  are  shown  in  Figures  44  and  45 .  In  both  cases, 

i 

thd  trend  is  found  to  be  toward  a  reduction  of  within-lot  variance.  This 

i 

indicates  the  overall  improvement  in  production  and  testing  processes.  Twin 

;  j 

City  shows  a  larger  reduction  as  compared  to  Lake  City. 

i  .  1 

•  I 

1  ’  •  , 

i  ! 
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CUMULATIVE  SUM 
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4 ,  CUTOFF  DATE  AND  CONTROL  LIMITS 

4 . 1  Method  ox  Approach 

Present  control  limits  for  ammunition  production  are  wide  leading 
to  severe  gun  design  requirements.  Attempts  are  therefore  to  be  made  to 
narrow  the  control  limits . 

It  is  observed  that  the  chamber  pressure  plots  show  a  downward 
trend.  This  can  be  attributed  to  drastic  process  changes  or  the  initial 
experimental  stage.  If  the  chamber  pressure  control  limits  are  based  upon 
chamber  pressure  values  after  the  process  changes  are  completed  (i.e., 
after  the  cutoff  date) ,  narrower  control  limits  can  be  obtained. 

Cumulative  sum  technique  vdth  reference  constant  equal  to  the 
mean  of  the  series  is  used  to  determine  the  cutoff  date.  An  empirical 
approach  using  cumulative  distribution  function  is  then  employed  to  obtain 
the  control  limits.  These  control  limits  are  based  upon  data  after  the  cut¬ 
off  date . 

4 . 2  Cutoff  Date 

The  cutoff  date  should  satisfy  the  following  requirements: 

(1)  It  should  be  possible  to  narrow  the  control  limits  when  based  upon 
observations  after  the  cutoff  date.  This  implies  that  the  chamber 
pressure  data  after  the  cutoff  date  should  be  sufficiently  smaller  in 
magnitude  than  the  data  before  the  cutoff  date . 


(2)  This  reduction  in  magnitude  should  be  maintained  for  a  considerable 
past  period  of  time  so  that  it  could  be  attributed  to  an  improvement 
in  process  rather  than  to  cyclic  or  chance  variations. 

(3)  Information  might  be  available  regarding  the  date  on  which 
process  changes  intended  to. reduce  the  chamber  pressure  were 
introduced.  If  the  observed  point  of  shift  corresponds  to 
this  prior  information,  then  it  may  be  taken  as  the  cutoff 
date. 

Method  employed  to  determine  the  cutoff  date  is  as  follows:  Com¬ 
puter  plots  of  available  chamber  pressure  data  are  obtained.  Th*  plots  are 
updated  as  more  data  become  available.  Usually  the  visual  picture 
indicates  whether  a  point  of  shift  exists;  for  example,  Figures  46  and  47 
show  the  tracer  chamber  pressure  plots  for  Lake  City  and  Twin  City  respec¬ 
tively.  (These  plot&  are  updated  over  those  given  in  Figures  20  and  21). 

A  careful  examination  of  these  plots  shows  a  point  of  shift  approximately 
80  lots  away  from  the  initial  reference  value  in  both  cases.  No  such  down¬ 
ward  shift,  is  observed  in  the  data  from  Winchester  (Fig.  48). 

Cusum  plotting  is  now  employed  to  quantify  the  point  of  shift.  The 
resulting  cusum  plots  for  Ball  and  Tracer  ammunition  from  different  manu¬ 
facturers  are  shown  in  Figures  49  through  56.  It  is  assumed  that  these 
points  of  shift  satisfy  the  conditions  necessary  for  cutoff  date.  The  point 
ot  shift  gives  the  lot  number  where  the  shift  occurred.  Cutoff  date  is  the 
corresponding  date  given  in  the  standard  lot  testing  results.  Cutoff  dates  are 


tabulated  in  Table  4 . 
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Figure  47  - - - - — 

chamber  pressure 
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4 . 3  Control  Limits 

Control  limits  are  based  upon  data  after  the  cutoff  date.  The 
following  procedure  is  used  to  determine  control  limits . 

The  range  of  chamber  pressure  readings  (range  =  maximum 
chamber  pressure  -  minimum  chamber  pressure)  is  subdivided  into 
twenty  equal  parts  and  the  number  of  lots  belonging  to  each  sub¬ 
group  is  calculated .  These  numbers  are  then  divided  by  the  total 
number  Oi  lots  to  obtain  an  estimate  of  the  probability  of  belong¬ 
ing  to  each  subgroup.  Cumulative  sum  of  these  probabilities  is 
plotted  against  the  corresponding  chamber  pressure.  The  resulting 
plo*.  is  known  as  the  empirical  cumulative  distribution  function 
(empirical  c.d.f.) . 

Only  one  sided  (upper)  control  limit  need  to  be  calculated  for 
the  chamber  pressure  data.  This  is  taken  as  the  99  percentile  of 
the  empirical  c.d.f. 

Plots  of  empirical  c.d.f.  for  chamber  pressure  data  (Ball  and 
Tracer)  from  the  different  manufacturers  are  given  in  Figures  5  7 


through  64 .  The  calculated  control  limits  are  given  in  Table  5 . 
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TABLE  4 

Cutoff  Date  for  Ball  and  Tracer  Ammunition 


Manufacturer 

Cutoff  Date 

Ball 

Tracer 

Lake  City 

June  21,  1968 

y  -  'iuary  30,  1970 

Twin  City  % 

March  10,  1969 

irril  15,  1969 

Winchester 

February  22,  1967 

Nil 

Remington 

Nil 

y 

Federal 

Nil 
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TABLE  5 


Control  Limits  for  Ball  and  Tracer  Chamber  Pressure 


Manufacturer 


Control  Limits 


Ball 

Tracer 

Lake  City 

49700  psi 

50500  psi 

Twin  City 

49200  psi 

51200  psi 

Winchester 

51000  psi 

51600  psi 

Remington 

49800  psi 

X 

Federal 


50600  psi 


X 
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f-.",: ; .  ,  5  .  AMMUNITION  SE LEGTION  FORWEAPON,. TESTS ,  ■  • .  - ,  • 

* 

5 ,1  Method  of  Approach 

The  capabilities^, of  contemporary  ammunition  production  are  such 

that  differences  exist  between  different  ammunition  lots.  Gun  perfor¬ 
ms?  •  -v.v  ■  o’  *  ’  1  "  ■  "  1  . 

mance  is  influenced  by  ballistic  properties  of  ammunition,  therefore, 

-<>■  •  •  i't  ■  ' «  •  ■  -i  »  •.*  -  •  ■  •  -  .  •  .  ■ 

selection  of  proper  ammunition  lots  for  weapon  testing  is  of  importance  . 

The  ammunition  responses  that  determine  gun  performance  are  not 
yet  known.  The  present  analysis  assumes  that,  gun  performance  is  governed 
by  chamber  pressure  and  port  pressure.  The  analysis  of  ammunition  pro¬ 
duction  process  indicates  the  presence  of  a  cutoff  date.  Since  the  consider¬ 
ation  of  ammunition  lots  from  the  production  period  before  the  cutoff  date  leads 
to  more  severe  gun  requirements ,  ammunition  lots  from  the  production  period 
after  the  cutoff  date  alone  have  been  used  to  develop  the  selection  criterion. 

The  analytical  approach  is  to  use  »  probabilistic  viewpoint  in  select¬ 
ing  those  ammunition  lots  from  normal  production  that  are  likely  to  give  a 
large  variation  in  weapon  performance.  The  analysis  is  first  detailed  for 
ammunition  from  Lake  City  Ammunition  Plant  and  a  similar  analysis  is  then 
carried  out  for  Twin  City  Ammunition  Plant.  The  cutoff  date  for  Lake  City 
is  July  21st  1968  and  for  Twin  City  is  October  3rd  1969,  The  analysis  for 
Lake  City  is  based  upon  a  total  of  942  lots,  covering  a  period  from  July  21st 
1968  to  January  1st  1971.  A  total  of  411  lots  covering  a  period  from 
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October  3rd  1969 


from  Twin  City. 


.  fo  • ' 'i orjA  c  * , 


5 . 2  Selection  Ba&ed(1Upon  ‘Chamber  PreBBarei(I&ke  Cl 

-tolieq  rtiiyj  .  aaoi  i.oi  1'Hdii‘fqL,  -‘utfisl'i:  .ov,i-  .  jxtr  - 

Chamber  Pressure  from  standard  tests  is  assumed  to  control  the 

xtftoiorerii  .fic.r}rrt?rr  >-u  io  -  r**-. >*«<. ic  f  .■» 

gun  performance.  Figure  65  is  the  plot  of  Lake  City  chamber  pressure 

.  sV’^Aicami  io  8i  g'l:'*;  r»i  no’*;C-  •  .  -  -  •  a*'  1  ■ 

against  the  corresponding  ammunition  lots  covering  the  period  after  the 


■»  «.  i 


cutoff  date.  The  range  of  chamber  pressure  readings  (range  =  max. 

hs-i-'.  ^  Ji>  -L  \  '  •.  ■  .  ' 

chamber  pressure  -  min.  chamber  pressure)  is  subdivided  into  twenty 

5  •>..  .  :  ,  ' 

equal  parts  and  the  number  of  lots  belonging  to  each  subgroup  is  calcula- 

-  3»~  '  >  ;.,H  V,ji . ,  .4  '  ■<  r 

ted.  Figure  66  shows  the  resulting  histogram.  Figure  67  is  a  plot  of  em- 

•iti&jfl  %:?.■  i-V.  ■■  ,  .  *■  .  ■■‘3, 

pirical  cumulative  distribution  function  (empirical  c.d.f.).  This  is  ob- 

0  >  ‘  VI  H  lv*,  JO  9!  •  ■  ...  -i  • 

tained  by  first  calculating  the  probability  of  belonging  to  each  subgroup. 


Empirical  c.d.f.  is  a  plot  of  the  cumulative  sum  of  these  probabilities, 
plotted  against  the  corresponding  chamber  pressure.  The  entire  sequence 
of  calculations  is  tabulated  in  Table  6 . 

’«  chamber  pressure  controls  the  gun  performance,  then  a  large 

*,  t  3  it*  *  » •  1  .  . 

variation  in  gun  responses  would  be  obtained  by  selecting  lots  with  chamber 

t/ J*  *  >  3,1  k\  t’j  .  v 

pressure  toward  the  high  and  the  low  ends  of  the  spectrum.  To  these  cate- 

gories  of  High  and  Low  lots  another  category  of  Medium  lots  is  added  to 

■•a '  s  y'  >  n  -  ■  i>‘*  •  ■  <  '  •  «  ■ 

determine  gun  performance  that  will  be  most  frequent  ly  met  in  practice! . 


Figure  65 _ 
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HISTOGRAM  AND  CUMULATIVE  DISTRIBUTION  FUNCTION 
LAKE  CITY  CHAMBER  PRESSURE 


Serial  ,  Number  of  Cumulative 

Number  Interval  Occurences  Probability  Probability 

\  _ 

1 


1 

372.00 

to 

378.55 

1. 

.001 

.001 

2 

378.55 

to 

385.10 

0. 

.000 

.001 

3 

385.10 

to 

391.65 

1. 

.001 

.002 

4 

391.65 

to 

398.20 

0. 

.000 

.002 

5 

398.20 

to 

404.75 

0. 

.000 

.002 

6 

404.75 

to 

411.30 

1. 

.001 

.003 

7 

411.30 

to 

417.85 

0. 

.000 

.003 

8 

417.85 

to 

424.40 

1. 

.001 

.004 

9 

424.40 

to 

430.95 

0. 

.000 

.004 

10 

430.95 

to 

437.50 

9. 

.009 

.014 

11 

437.50 

to 

444.05 

48. 

.051 

.065 

12 

444.05 

to 

450.60 

97. 

.103 

.168 

13 

450.60 

to 

457.15 

204. 

.216 

.384 

14 

457.15 

to 

463.70 

186. 

.197 

.582 

15 

463.70 

to 

470.25 

177. 

.188 

.769 

16 

470.25 

to 

476.80 

105. 

.111 

.881 

17 

476.80 

to 

483.35 

69. 

.073 

.954 

18 

483.35 

to 

489.90 

28. 

.030 

.984 

19 

489.90 

to 

496.45 

9. 

.009 

.994 

20 

496.45 

to 

503.00 

6. 

.006 

1.0 

Total  =  1 


Total=942 


The  classification  of  lots  into  High,  Medium,  and  Low  categories 


is  done^acpording  jtOf  the.  fpUpwijtigs  probabilistic  poitjt  of  view;.  .A  lot  is 
cla  s.sifiedc.ascHigh  .ifthe}  probabijityjof  getting  a,  chamber  ,pres  surer  higher 
than  the  chamber  , pressure  . of  that  particular  lot  is.  5% .  (.Similarly,,  if  the 
probability  of?  getting}  a  chamber  pressure  lower  thanjthat  of  a  particular 
lot  is  5%,  the  lot  is  classified  as  Low.  Medium  lot  is  one  that  is  most 
probable.  High  and  Low  lots,  therefore,  correspond  to  the  95  and  5  per¬ 
centiles  of  the  empirical  c.d.f.  Medium  lot  corresponds  to  the  mode  of 
the  histogram.  Another  statistic  that  could  be  used  to  classify  Medium 
lots  is  mean. >  However,  mode. appears  more  appealing,,  because  it  would 
lead  to  the  determination  of  most  frequent  gun  performance  rather  than  the 
average  gun  performance  as  determined  by  the  mean. 

The  calculated  critical  values  of  chamber  pressure  are  given  below: 

High:  48300  psi 

Medium:  45488  psi 

Low:  44150  psi 

Table  7  gives  the  selected  ammunition  lots.  These  lots  have  chamber  pressure 
values  close  to  the  critical  values. 

5.3  Selection  Based  Upon  Port  Pressure  (Lake  City) 

The  lot  selection  can  also  be  based  solely  on  port  pressure,  where 
it  is  assumed  that  port  pressure  alone,  governs  weapon  performance. 
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Figure  68  'Is  the'  plot  of ‘port  pressure  against  Ammunition  lot&;'  covering 
the1  period7  of  ammunition  production  afteMhb  cutoff  date.  Figure  69  is 
the  his'togram  of  port  pressure  and  -Fig-.-  70  is  the  empirical  c.d.f. .  •  The 
calculationsJare  detailed  In  Table  8.  '  The  critical  values  of  port  pressure 
for  classification  into  the  three- groups5 are  as  follows: 

High:  15921  psi 

Medium:  15430  psi 

Low:  14738  psi 

Table  9  gives  the  selected  lots  that  have  port  pressure  values  close  to  the 
critical  values. 

5 .4  Discussion  of  the  Two  Methods  of  Ammunition  Selection 

Both  the  methods  of  selection  are  extremely  simple  to  use,  however, 
they  suffer  from  the  following  disadvantages: 

(1)  A  comparison  of  Table  7  and  Table  9  reveals  that  classification 
based  upon  chamber  pressure  alone  and  based  upon  port  pressure 
alone  leads  to  different  lot  selection.  For  example,  lots  L- 1-2 04 
and  L-l-158  are  classified  as  Medium  based  upon  the  port  pressure 
but  are  classified  as  Medium  and  Low  respectively,  based  upon  the 
chamber  pressure.  It  is,  therefore,  necessary  to  determine  whether 
it  is  the  chamber  pressure  or  it  is  the  port  pressure  that  controls  gun 


performance . 


.'/raVy i  *3  nr  ? >" ‘jAiJr? 4  &  rS’  *il*>r*v ' 


T&BLE  7 


(Based  on  Chamber  Pressure  Alone) 
Lake  City 


Medium  Low 


Lot 

No. 

Chamber 

Pressure 

psa. 

Lot 

No. 

Chamber 

Pressure 

psi 

L- 

L- 

1-288 

45400 

1-284 

44200 

1-287 

45300 

1-79 

44200 

1-275 

45400 

1-177 

44100 

1-216 

45300 

1-174 

44100 

1-211 

45400 

1-158 

44200 

1-204 

45300 

1-1,42 

44300 

8SURE 


SQC.O  7-0Q...0.  -SCC.C  900. 


NUMBER  OF  OCCURENCE'S 
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Figure  :  69 
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TABLE  8 

HISTOGRAM  AND  CUMULATIVE;  DISTRIBUTION  FUNCTION 
LAKE  CITY  PORT  .PRESSURE  . 

’  _/2‘-  J-  r  ;•  t:rrr? rm  ■■• 111  ■•  r- rr — r-rr-rrer-** -**• 


Serial  Number  of  Cumulative 

Number  Interval  Occurences  Probability  Probability 


1 

1300.00 

to 

1318.00 

1. 

.001 

.001 

2 

1318.00 

to 

1336.00  ' 

0. 

.000 

.001 

3 

1336.00 

to 

1354.00 

0. 

.000  \ 

.001  ’ 

4 

1354.00 

to 

137.2.00 

2/. 

.002 

.003, 

5 

1372.00 

to 

1390.00 

3. 

.003 

.006' 

6 

1390,00 

to 

1408.00 

2. 

.002 

.008' 

7 

1408,03’ 

to 

1426.00 

3. 

.003 

.012 

8 

1426.00 

to 

1444.00 

8. 

•  .008 

.020 

9 

1444.00 

to 

1462.00 

13. 

.014 

.034 

10 

1462.00 

to 

1480.00 

35. 

,037 

.071 

11 

1480.00 

to 

1498.00 

41. 

.044 

.115 

12 

1493.00 

to 

1516.00 

146. 

.155 

.270 

13 

1516.00 

to 

1534.00 

217. 

.230 

.500 

14 

1534.00 

to 

1552.00 

236. 

.250 

.750 

15 

1552.00 

to 

1570.00 

135. 

.143 

.894 

16 

1570.00 

to 

1588.00 

44. 

.047 

.940 

17 

1588.00 

to 

1606.00 

41. 

.043 

.984 

18 

1606.00 

to 

1624.00 

12. 

.013 

.997 

19 

1624.00 

to 

1642.00 

2. 

.002 

.999 

20 

1642.00 

to 

1660.00 

1. 

.001 

1.0 

Total»942  Total=l 
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Lot  Selectxon—yBasea  on-Port* Pressure  Alone) 


_ _  _ Lake  City.  . 

.■  fjj-dimuC  p 
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Port  >• 
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'  Port 
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Pressure;  •, 

Lot  ■ 

Pressure, 

r  ,  LOt  - 

Pressur 

No. 

soft.... 

psi 

T  t\f  i 

No.  ‘ 

psi 

No. 

„  .  isi 

l-300  , 

*  O') . 

L- 

L-, 

v  , 

1-3,18 

15870  j . r, 

1-359 

15410 

1-297 

.14800 

1-317 

15970  LW 

1-338 

15450  ' 

1-287 

14900 

1-308 

'* f  0 . 

•  'OC 

15940 

’*  l  0 

1-296 

1539b' 

1-2  8  5 

"  14720 

*«  u  ; 

1729,8 

15970 

1-281 

15410, 

1-261 

,14870 

1-247 

15930  ■' 

1-274 

15440 

:  1-217 

*  14830 

*  '*  4 

cc. 

• 

1-265 

15410 

i, 

< '  •  i 

•  . 

1-250 

15420. 

i 

{  :  , 

; 

1-239 

15430 

• 

V  *'  <v  , 

»i  ^  , 

1-204 

15400 

»  * 

fu, 

1-158 

1543,0 

109 


* n  r 


(2)vWGvm,j>e.rfqnnance,  jls>  more  4ikalyatof.be  controlleduby>the  parameters  ; 
r  .6qf,<pia|ip  pressurertime  curv&  :It;  is  ipoasibiethat,  for>different-  . 

...  ^prassurextime  curyesvthei.sama  chamber,  pres  sure  rvalue  4s  obtained 
Sameftcan:.ibe  »spid  (regarding  ithei  port  pressure .  It  ;ie p  therefore;  ^rather 
•,r  unlikely  ;that,gun  <parfprmancewould'  beidetermined^by chamber  pressure 
alone, (or,, port  prqssurp}  alone/; „ ;  ,  •;  c  \-U 
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5.5  Selection  Based  Upon  Chamber  Pressure  and  Port  Pressure  (Lake  City) 

■,  ’  ■  ,u<[  ’•  *•  i'C  '  -  >  --'o'  *  •• 

It  is  considered  more  probable  that  gun  performance  is  governed 

jointly  by  the  chamber  pressure  and  the  port  pressure.  The  final  selection 

.  f:*/v,r  i  y  '  *'  ‘  n'‘ 

is.  therefore,  based  upon  both  the  chamber  pressure  and  the  port  pressure. 


Table  104s.  the  bivariate  histogram  of  chamber  pressure  and  port  pressure 


excluding  the  outlays.  This,  computer  generated  table  gives  the  number  of 
lots  belonging  to  a  certain  chamber  pressure  range  and  a  certain  port  pressure 
range.  For  example,  there  are  26  lots  that  have  a  chamber  pressure  between 
45.800  psi  and  46100  psi  and  a  port  pressure  between  15400'psl  and  *15500  psi. 
In  this  table  the -group  width  (  A^)  along  the. port  pressure  axis  is  equal 
to  half  standard  deviation  (S.D.)  and  alongthe  chamber  pressure  axis 
(  Ag)  is  equal  tc  one  S.D-.  , of  chamber ipressure .  (The  mean  value  of  S.D. 
is  1330  psi ^ence,. S.D, « ot mean  chamber  pressure  =  1330/  20 a  300  psi). 

The  selected  frequencies  are  enclosed  in  small  rectangles  in  Table  10. 


Medium  lots  are  the  most  probable,  hence  the  central  rectangle  is  so  placed 


no 


m 


as  taxrOvei^the. makimui^HufriberofilotssfeTfe^ 

lots  is  prftnnar3a^9oMrobdo&^4ll«^hflttiitber^[^^  ThOiTe^Kave 

c(2  f^rSa  D1.  )>  as  *thet  Medium'  lots'  bat^tli#  <sham- 
l&gipteest&foradge  dfe  ^situatedearauhd  <the  5  high  <andlow'  criticaFvalues 
^detemtaadfbyrOsingichamber^pressiure'crlterlbWaloWa  i-  Tor  'iPlonhi^  Selec¬ 
tion'  the  chamber  pressure  range  is  2  S .  Dv^and-for'-Tight  Selection'’ 


by  shaded  rectangles)  the  chamber  pressure  range  is  1  S.D. . 

(',*74  ’«  i)  -Yil?  *??'.  .’f’.  *■£  .  ' l{ .  ‘  'y ,  '  t 

Table  X  gives  the  critical  values  of  chamber  pressure  and  port  pressure  in 


'  t  *  '  '  TT '  *'  C  s'  . 

both  the' cases.  Table  12  gives  the  selected  lots  for  Tight  Selection  and  the 
selected’ lots' ''for  ‘Normal' 'Selection  are  given  in  Table  13. 


.*«.  >w<v«  »  i.t  }m  aw.  a  <>>“?  ^ 
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5 ,6  Selection  Based  Upon  Chamber  Pressure,  and?  Port  Pressure  (Twin  City) 


rsr 


7o  -  s\-I  <~/L  similar  approach  isnowused-forTwin  Gity  datav^The  Twih  Cify 
«  acharaber  pre  s  8ure;andriport;pre  s  sure'for;  the  period  of*  normal  production  are 
»fiploKed  inyFigures^‘71 4shd  72. *  Fig;  73  gives  the  histogram  of  Twin  City 
j^chaadbOripreSsure^nd  Fi'jv*74;4a  the  plot1  of  empirical  c.d'.f. .  The  entire  ' 
sequence -of  icalculations;  is  shown  in  Table  14 .  The  bivariate  histogram 
of  chamberapressure^and  ‘port  pressure  is 'shown  in  Table  15.  Table  16  in¬ 
dicates'  thaicriticali values  for  Normal  and  Tight  Selection.  The  selected' 

.  tilot^for  Tight  Selection  and  Normal  Selection  are  included  in  Tables  17  and 
.uilSirejSpoctivOly^  ^ 
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BIVARIATE  HISTOGRAM  OF  CHAMBER  PRESSURE  AND  PORT  PRESSURE 
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Note:  Values  of  pressures  along  the  two  axes  are 
maximum  values  for  that  group. 
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■TABLE  14 

HISTOGRAM  AND  CUMULATIVE  DISTRIBUTION  FUNCTION 

twiiTcity  cham&er  pressure  ' 


Serial  Number  of  Cumulative 

Number  Interval  Occurences  Probability  Probability 


1 

415.00 

to 

419.50 

2. 

.005 

.005 

2 

419.50 

to 

424.00 

2. 

.005 

.010 

3 

424.00 

to 

428.50 

1. 

.002 

.012 

4 

428.50 

to 

433.00 

9. 

.022 

.034 

5 

433.00 

to 

437.50 

6. 

.015 

.049 

6 

437.50 

to 

442.00 

19. 

.046 

.095 

7 

442.00 

to 

446.50 

17. 

.041 

.136 

8 

446.50 

to 

451.00 

44. 

.107 

.243 

9 

451.00 

to 

455.50 

34. 

.083 

.326 

10 

455.50 

to 

460.00 

50. 

.122 

.448 

11 

460.00 

to 

464.50 

65. 

.158 

.606 

12 

464.50 

to 

469.00 

42. 

.102 

.708 

13 

469.00 

to 

473.50 

40. 

.097 

.805 

14 

473.50 

to 

478.00 

31. 

.075 

.881 

15 

478.00 

to 

t  50 

12. 

,029 

.910 

16 

482.50 

to 

4b/. 00 

17. 

.041 

.951 

17 

487.00 

to 

491.50 

10. 

.024 

.976 

18 

491.50 

to 

496.00 

6. 

.015 

.990 

19 

496.00 

to 

500,50 

3. 

.007 

.998 

20 

500.50 

to 

505.00 

1. 

.002 

1.0 

Total=411  Total=l 
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1160  15140 


(Based  Upon  Chamber  Pressure  and  Port  {Pressure) 
Chamber  Pressure  Tolerance  =  2  S.D. 
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5 .7  Design  of  Experiments  to  Obtain  the  Effect  of  Ammunition 
Characteristics  on  Weapon  Performance 

In  selecting  ammunition  for  weapon  tests,  it  has  been  assumed 

that  the  weapon  performance  is  governed  by  the  copper  crusher,  chamber 

j  c  r>  o  .'•>  -  ^  <  . 

pressure  and  port  pressure  characteristics  of  the  Ammunition*  lot..  This 

;  t'f  %  <  .■  ’*  -'i  ,  ■■  ->  | 

,  %  r-  >  *  ’  '  »  ■  •  •"  t  n.  ^  * 

is  not  known  to  be  true.  ;The  true  ammunition  parameters  that  control, 

)  *  ^  «*'  l  *  j 

gun  performance  need  to  be  determined .  ^ Once.' these  parameters  are* known, 

i  '• 

•  \  I  M 

they  can  be  substituted  for  chamber  pressure  and. port  pressure;  and  the 

.  |  tf**1’*  **  *  > 

analysis  repeated  to  obtain  prqper;;ammunitipn*'lot8  for  weapon  tests. 


Lacking  a  complete  theoretical  analysis  at  the  present  time,  true 

I  «  >  « 

ammunition  parameters  can  only  be  determined  experimentally,.  To  minimize. 

t 

experimentation,  proper  experimental  design  techniques  have  to  be  used. 
One  such  experimental  procedure  is  now  illustrated . 


*  ^ 

Let  us  suppose  that  the  effect  of  chamber-pressure, and  port  pressure 

{  2 

on  cyclic  rate  is  to  be  determined .  Figure  75  shows  a -replicated  2.  factorial 
design  with  cyclic  rate  as  response  .  For  tests  at  each  of  the  four,  points,  two 
ammunition  lots  are  selected  with  their  chamber  pressure  and  port  pressure 

j,- 

within  the  indicated  ranges.  Since  these  ranges  are  twp  standard  deviations 


wide,  the  lots  are  considered  identical.  The  two  tests  at  each  of  the  four 


points  are,  therefore,  replicates.  This  design  enables  the  determination  of 


the  main  effects  of  chamber  pressure  and  port  pressure  as  well  as  the  effect 
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Such  factorial  and/or  fractional  factorial  designs  can  be  employed 
to  determine  ammunition  parameters  that  show  significant  effects  on  weapon 
responses .  These  parameters  can  then  be  said  to  control  weapon  perfor¬ 
mance  and  would  be  used  to  select  ammunition  for  weapon  testing. 
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6.  DETERMINATION  OF  PROPER  NUMBER  OF  TESTS 
FOR  ESTIMATING  AMMUNITION  PARAMETERS 

6 . 1  Method  of  Approach 

To  obtain  good  estimates  of  ammunition  parameters  such  as 
chamber  pressure,  port  pressure,  peak  chamber  pressure,  velocity,  etc., 
proper  number  of  tests  have  to  be  conducted.  If  smaller  than  necessary 
number  of  tests  are  conducted,  the  parameters  would  be  estimated  with 
large  variance .  On  the  other  hand ,  conducting  more  than  necessary 
number  of  tests  involves  unnecessary  expenditure.  There  is  a  need  to  ob¬ 
tain  an  appropriate  number  of  tests  to  be  performed . 

This  number  would  be  determined  by  the  desired  precision  with  which 
the  parameters  are  to  be  estimated  and  the  experimental  error  involved.  For 
example,  in  the  measurement  of  chamber  pressure  by  standard  tests  if  the 
desired  precision  is  r  =  1  and  the  experimental  error  is  known  to 

V 

be  o2,  then  the  required  number  of  tests  'n*  to  get  tho  desired  precision 
is  given  by 

2,2 

n  »»  o  /  Oj 

The  chamber  pressure  is  then  estimated  by  the  mean  of  these  n 

2 

readings.  The  moan  has  an  error  variance  of  which  is  a  prespacified 
value.  The  only  unknown  quantity  is  o2  —the  error  variance  for  each  test. 
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Hence  a  has  to  W  expert merits lly  determined'.  The  experimental  error  is 

r*  f  *  '  *•  V  .*  -  •  >• 

•*  *  ‘  t  ■  A 

a  function  of  the  foliowing: 

(1)  The  parameter  or  parameters  to  be  estimated,  i.e.,  the  error  in¬ 
volved  in  the  measurement  of  chamber  pressure  can  be  different 
from  the -error  in-the  measurement  of  port  pressure. 

(2)  Experimental  setup  which  includes  variables  like  barrels,  measur¬ 
ing, devices  (copper  crusher,  piezo)  cartridges  from  different 
manufacturers  etc. 

(3)  Inherent  variability  of  the  manufacturing  process  in  producing  a  lot 
of  ammunition.  Table  Vindicates  that  lots  belong  to  the  HIGH 

,r  .category  have  a  somewhat  larger  chamber  pressure  standard 
deviation  as  compared  to  the  lots  belonging  to  the  LOW  category. 

If  the  population  of  ammunition  lots  is  considered  to  be  homogeneous, 
then  the  average  value  of  standard  deviation  of  chamber  pressure  (e.g.  1330 

psi  for  Lake  City)  is  a  good  estimate  of  standard  error  in  copper  crusher 

-  / 

testing,  for  the  measurement  of  chamber  pressure.  However,  differences 
exist  amongst  the  ammunition  lots.  It  is,  therefore,  necessary  to  conduct 
experiments  to  determine  the  experimental  error. 

6.2  Sequential  Procedure  for  Estimation  of  Experimental  Error 

A  sequential  procedure  is  now  illustrated  by  which  an  estimate  of 
experimental  error  can  be  obtained.  Figures  76,  77,  and  78  are  based  upon 
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20  individual  copper  crusher  chamber  pressure  readings  per  lot.  Tho  dote 

*  *  'Y 

were  obtained  from  B.A.A.P.  Figures  79  end  80  are  based  upon  10  individ¬ 
ual  piezo  peak  pressure  readings  obtained  from  Frankford  Arsenal.  Tho  en¬ 
tire  set  of  data  is  given  in  Appendix  III. 

The  figures  illustrate  sequential  estimation  of  variance  of  obser¬ 
vations  using  2,3,4,  ...  number  of  tests.  As  can  be  observed,  the  estimate 
of  variance  fluctuates  considerably  in  the  beginning  but  tends  to  stabilize 
as  the  number  of  observations  increases.  The  value  about  which  the  variance 
stabilises  is  the  estimate  of  experimental  error.  The  figures  indicate  that  in 
all  coses  the  variance  has  not  stabilised .  The  conclusion  is  that  20  readings 
are  not  sufficient  to  estimate  experimental  error  in  copper  crusher  testing 

and  same  is  true  about  ten  readings  in  piezo  testing. 

It  should  be  noted  that  the  current  practice  in  standard 
acceptance  testing  is  to  conduct  20  tests  to  obtain  an  estimate 
of  standard  deviation.  In  view  of  the  analyses  conducted  above, 
this  estimate  of  standard  deviation  has  a  large  variance  asso¬ 
ciated  with  it.  It  is  suggested  that  sufficiently  large  number 
of  tests  be  conducted  to  obtain  a  stable  estimate  of  experimen¬ 
tal  error. 

Once  such  an  estimate  of  experimental  error  is  available, 
necessary  number  of  tests  to  obtain  parameter  estimates  with 
desired  degree  of  precision  can  be  determined. 
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7.  ANALYSIS  OF  DATA  AT  DIFFERENT  STAGES  OF  MANUFACTURE 

Data  are  generated  at  different  stages  In  the  manufacturing  process . 

An  analysis  of  these  data  would  Indicate  changes  in  ammunition  character¬ 
istics  during  the  process  of  manufacture.  Such  an  analysis  would  help 
explain  the  final  ballistic  characteristics  attained  by  a  given  ammunition. 

7.1  Time  Series  Analysis  of  Propellant  Lot  Characteristics  from  B.A.A.P. 

Since  the  propellant  lots  are  produced  sequentially,  a  relationship 
may  be  expected  between  successive  lots.  The  ballistic  results  are  ob¬ 
tained  by  adjusting  the  charge  weight  to  obtain  muzzle  velocity  close  to 
3250  f.p.s.  Changes  in  lot  characteristics  would  be  reflected  in  the  set 
charge  weight  and  the  resulting  chamber  pressure.  These  two  responses, 
therefore,  have  been  used  to  represent  the  properties  of  the  propellant  lots. 

The  series  of  charge  weight  and  chamber  pressure  are  given  in  Appendix  IV. 

Plot  of  chamber  pressor®  data  is  shown  in  Fig.  81.  Fig.  82  shows 
the  plot  of  the  autocorrelation  function  and  the  partial  autocorrelation  function 
for  the  original  series.  A  plot  of  the  first  difference  of  the  original  series 
is  shown  in  Fig.  83  and  the  corresponding  auto  and  partial  correlation  functions 
are  shown  In  Fig .  84 .  The  plots  indicate  that  the  first  difference  of  the 
chamber  pressure  data  is  stationery  and  follows  an  MA  (1)  process.  The 
mode!  is,  therefore,  identified  asr 
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(1  -  0B)at  =  at  -  0 


•*  '  *  .  >  *  . 

where 

yt  ®  chamber  pressure  for  the  t*h  lot 
i  u.  ..  ,at.  *  NID  (0.O2) 

0  ®  parameter  to  &e  estimated 

The  value  of  0  Is  found  by  nonlinear  estimation  and  Is  0.8892, 
Fig,  85  shows  the  plot  of  the  residual  aj.  s.  The  model  is  found  to  be 
adequate. 

Fig,  86  through  90  show  a  similar  analysis  for  the  charge  weight 
series  and  the  time  series  model  for  this  data  is  found  to  be 


yt  ’  yt-1  =  (1  -  0.5425B)at 

where 

yt  =  charge  weight  for  the  tth  lot 
at~NID  (0,  o2). 


These  models  are  subject  to  the  following  limitationfr 

(1)  Testing  is  done  using  components  from  five  different  manufacturers. 
The  differences  in  components  are  likely  to  influence  the  results. 

(2)  Lots  do  not  strictly  represent  the  time  sequence  of  powder  production 
process  because  they  are  obtained  by  blending  fifteen  preblends  not 
necessarily  in  the  order  in  which  the  preblends  were  produced. 
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Within  these  limitations,  the  successive  propellant  lots  can  be 
said  to  exhibit  dependant  characteristics.  If  it  can  be  assumed  that  the 
successive  propellant  lots  are  used  at  the  ammunition  plants,  the  depen- 
dencybe  tween  propeilant  lots  would  cause,  to  a  small  extent,  the  suc¬ 
cessive  ammunition  lots  to  show  correlated  characteristics. 

7.2  Comparison  of  Acceptance  Test  Results,  for  the  Same  Propellant  Lots, 
from  B.A.A.P. ,  Lake  City  and  Twin  City 

Each  propellant  lot  produced  at  B.A.A.P.  is  in  effect  tested  for 
acceptance  twice,  once  at  B.A.A.P.  and  then  at  the  place  of  cartridge 
manufacture .  There  is  a  time  lag  of  a  few  months  between  the  two  testing 
dates  -  If  it  could  be  assumed  that  the  testing  equipments  are  alike,  then  a 
comparison  of  the  two  results  would  indicate  the  changes  in  powder  proper¬ 
ties  during  the  elapsed  time. 

Table  19  gives  the  propellant  lot  acceptance  test  results  from 
B.A.A.P. ,  conducted  with  Lake  City  components.  Table  20  gives  the  cor¬ 
responding  results  from  Lake  City.  Table  21  is  the  subtraction  of  Table  20 
from  Table  19.  Tables  22,  23  and  24  give  similar  results  with  Twin  City 
components.  It  can  be  observed  that  the  quantities  in  Tables  21  and  24  have 
small  magnitudes  and  random  variations  in  signs.  The  conclusion  is  that 
powder  properties,  in  terms  of  the  responses  considered,  do  not  seem  to 
change  during  the  period  of  few  months  between  the  two  testing  dates. 
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Figure  89 (a) 


TABLE  19 


PROPELLANT  LOT  ACCEPTANCE  TEST  RESULTS  FROM  BAAP 


Date 

Fired 

Propellant 
Lot  No. 

Charge 

Weight 

(Grains) 

Corrected 

Velocity 

(fps) 

Corrected 

Chamber 

Pressure 

(psi) 

Corrected 

Port 

Pressure 

(psi) 

8-14-69 

46362 

27.9 

3244 

45.100 

15500 

8-18-69 

46364 

27.9 

3243 

46700 

15700 

8-19-69 

46366 

27.9 

3240 

46800 

15400 

8-21-69 

46412 

27.7 

3255 

47500 

15500 

9-  3-69 

46414 

27.8 

3251 

46100 

15500 

9-10-69 

46416 

28.0 

3241 

46600 

15500 

9-15-69 

46418 

27.5 

3242 

44200 

15500 

9-22-69 

46425 

27.6 

3248 

44500 

15600 

9-30-69 

46427 

27.9 

3241 

45200 

15800 

10-10-69 

46430 

28.3 

3245 

47700 

15700 

10-28-69 

46881 

28.1 

3249 

48000 

16100 

10-30-69 

46883 

27.9 

3242 

45800 

16800 

11-  7-69 

46886 

27.6 

3242 

45900 

15200 

11-21-69 

46890 

27.1 

3240 

44700 

15700 

12-16-69 

46893 

26.8 

3250 

45200 

15600 

1-13-70 

46900 

26.9 

3245 

45700 

15200 

1-29-70 

46507 

27.5 

3252 

46900 

15300 

2-11-70 

46508 

27.3 

3244 

45300 

15600 

2-24-70 

46514 

27.6 

3248 

46400 

15700 

2-25-70 

46515 

27.2 

3247 

45600 

15700 

2-27-70 

46516 

27.0 

3247 

48800 

15100 

3-  9-70 

46517 

27.2 

3244 

47900 

15600 

TABLE  IS 


Cont'd 


Date 

Fired 

Propellant 
Lot  No. 

Charge 

Weight 

(Grains) 

Corrected 

Velocity 

(fos) 

Corrected 

Chamber 

Pressure 

(psi) 

Corrected 

Port 

Pressure 

(psi) 

3-23-70 

46520 

27.2 

3255 

47000 

15100 

5-12-70 

46529 

27.8 

3248 

45400 

15700 

6-17-70 

46607 

28.0 

3244 

46700 

15500 

6-25-70 

46609 

27.9 

3245 

45800 

15800 

7-  8t70 

46611 

27.4 

3255 

47000 

15500 

27.59 

3246 

46240 

15588 

TABLE  20 


154-B 


PROPELLANT  LOT  ACCEPTANCE  TEST  RESULTS 
FROM  LAKE  CITY  AMMUNITION  PLANT 


Date 

Fired 

1970 

Propellant 
Lot  No. 

Charge 

Weight 

(Grains) 

Corrected 

Velocity 

(fps) 

Corrected 

Chamber 

Pressure 

(psi) 

Corrected 

Port 

Pressure 

(psi) 

2-  6  to 
2-13 

46362 

27.9 

3237 

45440 

15600 

1-30-70 

46364 

27.95 

3252 

47000 

15400 

1-26  to 
1-30 

46366 

28.0 

3271 

47200 

15700 

3- 31  to 

4- 8 

46412 

27.8 

3244 

46022 

15500 

2- 26  to 

3- 3 

46414 

27,75 

3249 

45677 

15380 

2-13  to 
2-17 

46416 

27.9 

3223 

44700 

15310 

4-8  to 
4-15 

46418 

27.7 

3238 

44955 

15510 

4- 29  to 

5- 16 

46425 

27.8 

3242 

45444 

15300 

2-17  to 
2-24 

46427 

27.8 

3243 

44477 

15530 

3-5  to 
3-10 

46430 

28.1 

3243 

46360 

15350 

3-10  to 
3-13 

46881 

27.85 

3240 

45833 

15370 

5-7  to 
5-13 

46883 

27.8 

3257 

45625 

15440 

3-16  to 
3-24 

46886 

27.8 

4356 

45622 

15610 

3-24  to 
3-30 

46890 

27.0 

3240 

45211 

15320 

4-15  to 
4-21 

46893 

26.9 

3228 

45325 

15370 

4-22  to 
4-28 

46900 

26.9 

3250 

45566 

15440 

154-C 
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!  TABLE  20.  Cont *d 


Date 

Fired 

1970 

Propellant 
Lot  No. 

Charge 

Weight 

{Grains) 

Corrected 

Velocity 

(fps) 

Corrected 

Chamber 

Pressure 

(psi) 

Corrected 

Por 

Pressure 

(psi) 

5-14  to 
5-19 

46507 

27.5 

3234 

44488 

15446 

5-21  to 
5-26 

46508 

27.5 

3246 

45420 

15271 

5- 27  to 

6- 14 

46514 

27.7 

3245 

45400 

15387 

6-5  to 
6-11 

46515 

27.3 

3253 

45666 

15310 

6-12  to 
6-26 

46516 

27.1 

3237 

46033 

15116 

6-18  to 
6-24 

■  46517 

27.2 

3234 

46300 

15080 

6-29  to 
6-30 

46520 

27.1 

3250 

46636 

15150 

6- 30  to 

7- 13 

46529 

27.6 

3258 

46009 

15350 

7-14  to 
7-22 

46607 

27.9 

3248 

46212 

15300 

7-22  to 
7-27 

46609 

27.7 

3242 

46044 

15250 

7-27  to 
7-31 

46611 

27.3 

3243 

45588 

15480 

27.59 

3244.5 

45700 

15374 

DIFFERENCE  IN  BMP  RESULT  AND 
ME  Cl TI  AMMUNITION  PLANT  RESULTS 
(BMP  -  LAKE  CITY) 


Time 

Lag 

(Months) 

Propellant 
Lot  No. 

Difference 
in  charge 
weights 
(grains) 

Difference 

in 

Velocities 

(fps) 

Difference 
in  chamber 
Pressures 
(psi) 

Difference 
in  port 
Pressures 
(psi) 

6 

46362 

0 

7 

340 

100 

5  1/2 

46364 

-0.05 

-9 

-300 

300 

5  1/2 

46366 

-0.10 

-3?. 

-400 

-300 

7 

46412 

-0.10 

11 

1478 

0 

6 

46414 

+0.05 

2 

423 

120 

4 

46416 

0.1 

18 

1900 

190 

6 

46418 

-0.2 

4 

-755 

-10 

7 

46425 

-0.2 

6 

-944 

300 

4  1/2 

46427 

0.1 

-2 

723 

270 

5 

46430 

0.2 

2 

1340 

350 

4  1/2 

46881 

+0.25 

9 

2167 

730 

6  1/2 

46883 

0.1 

-15 

175 

1360 

4  1/2 

46886 

-0.2 

-14 

453 

-410 

4 

46890 

0.1 

0 

-511 

380 

4 

46693 

-0.1 

22 

-125 

230 

3  1/2 

46900 

0 

-5 

134 

-240 

3  1/2 

46507 

0 

18 

2412 

-416 

3  1/2 

46508 

-0.2 

-2 

-120 

329 

3 

46514 

-0.1 

3 

1100 

313 

3  1/2 

46515 

-0.1 

-6 

-66 

390 

3  1/2 

46516 

-0.1 

10 

2767 

-16 

3  1/2 

46517 

0 

10 

1600 

520 

3 

46520 

0.1 

5 

364 

-50 

2 

46524 

0.2 

-10 

-609 

350 

1 

46607 

-0.1 

-4 

488 

200 

1 

46609 

0.2 

3 

-244 

550 

1/2 

466H 

0.1 

12 

1412 

20 

I 
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TABLE  22 

PROPELLANT  LOT  ACCEPTANCE  TEST  RESULTS  FROM  BAAP 


Date 

Fired 

Lot 

Number 

Charge 

Weight 

(Grains) 

Velocity 

(fps) 

Corrected 

Chambor 

Pressure 

(psi). 

Corrected 

Port 

Pressure 

(psi) 

4-1-70 

46522 

27.4 

3257 

46100 

I 

15400 

4-22-70 

46524 

27.1 

3242 

46400 

15100 

5-4-70 

46527 

27.0 

3254 

45200 

15400 

5-7-70 

46528 

27.1 

3242 

48500 

1510Q 

5-22-70 

46530 

27.3 

3244 

44800 

13500 

1 

6-26-70 

46610 

27.7 

3245 

45500 

15800 

7-7-70 

46612 

27.4 

3254 

45200 

15700 

7-17-70 

46613 

27.2 

3245 

46600 

156Q0 

) 

7-28-70 

46617 

26.9 

3244 

46400 

15200 

7-2-70 

46618 

27.1 

3253 

45700 

15100 

8-4-70 

46620 

27.3 

3254 

46500 

15600' 

8-4-70 

46621 

27.8 

3249 

46000 

,  15900 

8-20-70 

46624 

27.6 

3247; 

45600 

15600' 

8-24-70 

46625 

27.5 

3247 

45500 

15700 

9-2-70 

46938 

27.4 

3252 

46800 

15600 

9-1-70 

46940 

27.5 

3  248 

44300 

15700 

9-17-70 

46942 

27.3 

3248 

44700  ' 

15200 

1  l 

9-23-70 

46944 

27.6 

3257 

46500 

15500 

9-30-70 

46945 

27.9 

3241 

45000 

15800  . 

9-30-70 

46946 

27.6 

3251 

47400 

15300 

1  -5-70 

46948 

27.4 

3250 

45200 

15400 

1 

10-9-70 

46949 

27.0 

3255 

43900 

15J300 

10-22-70 

46953 

27.1 

3240 

44500 

14900 

10-27-70 

46956 

27.4 

3254 

46100 

15600 

10-30-70 

46958 

27.2 

3251 

371*7*6 

46300 

TTTBTT 

15400 

134*6 

i 
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TABLE  2% 

PROPELLANT  LOT  ACCEPTANCE  TEST  RESULTS 
FROM  TWIN  CITY  AMMUNITION  PLANT 


i 


Date 

Fired 

1970 

Lot 

Number 

Charge 

Weight 

(Grains) 

Velocity 

(fps) 

Corrected 

Chamber 

Pressure 

(psi) 

Corrected 

Port 

Pressure 

(psi) 

6-1  to 
6-11 

46522 

27.4 

3243 

45920 

15484 

6-15  to 
6-22 

46524 

27.2 

3242 

45500 

15162 

6-22  to 
6-29 

46527 

27.3 

3246 

46180 

15432 

6- 30  to 

7- 9 

46528 

27.0 

3248 

46400 

15170 

7-12  to 
7-15 

46530 

27.4 

3244 

46100 

15207 

7-16  to 
7-27 

46610 

27.3 

3234 

45200 

15240 

7-30  to 
7-31 

46612 

27.3 

3252 

46150 

15360 

8-4  to 
8-12 

46613 

27.3 

3241 

45740 

15254 

8-13  to 
8-19 

46617 

27.1 

3232 

46475 

15150 

8-21  to 
8-31 

46618 

27.0 

3248 

45220 

15272 

8- 31  to 

9- 8 

46620 

27.3 

3241 

45325 

15375 

9-8  to 
9-16 

46621 

27.4 

3257 

46300 

15385 

9-17  to 
9-23 

46624 

27.8 

3237 

44425 

15648 

10-2  to 
10-8 

46625 

27.5 

3251 

45300 

15490 

9-24  to 
9-30 

46938 

27.4 

3244 

44575 

15815 

10-12  to 
10-15 

46940 

27.5 

3249 

44750 

15460 

156-B 


TABLE  23  Cont'd 


Date 

Fired 

1970 

Lot 

Number 

Charge 

Weight 

(Grains) 

Velocity 

(rps) 

Corrected 
Chamber ' 
Pressure 
(psi) 

Corrected 

Port 

Pressure 

(psi) 

10-16 

10-26 

to 

46942 

27.3 

3241 

46000 

15270 

10- 27 

11- 3 

to 

46944 

27.4 

3255 

46225 

15245 

11-4  to 
11-10 

46945 

27.7 

3248 

44400 

15605 

11-17 

11-18 

to 

46946 

27.7 

3248 

46220 

15382 

11- 19 

12- 1 

to 

46948 

27.5 

3243 

46400 

15205 

12-1  to 
12-11 

46949 

27.2 

3244 

45900 

15425 

12-12 

12-17 

to 

46953 

27.0 

3242 

46180 

15140 

12-21 

12-23 

to 

46956 

27.2 

3234 

46125 

15223 

12-28 

12-30 

to 

46958 

27.4 

3243 

47133 

15217 

27.34 

3244.28 

45765 

15342 

157 


i 
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TABLE  24 

i 

DIFFERENCE  IN  BAAP  RESULT  AND 
TWIN  CITY  AMMUNITION  PLANT  RESULTS 
(BAAP  -  TWIN  CITY), 


1  1 

Time 

Lag 

(Months) 

Lot 

Number 

Difference 
,  in  charge 
weight 
(grains) 

4 

Difference 

in 

Velocities 

,(fps) 

Difference 
in  corrected 
Chamber 
Pressure 
: (PSi) 

Difference 
in  corrected 
Port 
Pressure 
(psi) 

2 

46522 

,  0 

14 

180 

-84 

2 

46524 

-0.1 

1 

0 

900 

-62 

2 

46527 

,-0.3 

8 

-980 

-32 

2 

46528 

0.1 

-6 

2100 

-70 

1  1/2  , 

46530 

-0.1 

0 

-.1300 

293 

1 

46610 

,  0.4 

11 

300 

560 

1  i 

46612 

0.1 

2  , 

-950 

340 

1 

46613 

-0.1 

1  4 

860 

346 

1/2 

46617 

-0 . 2 

12  i 

-75 

50 

1 

466^8 

0.1 

'  5 

480 

-172 

1 

46620 

0 

13 

1175 

225 

1 

46621 

0.4 

-8  ; 

-300 

515 

1  i 

46624 

-0.2: 

10, 

1175 

-48 

2 

46625 

i  o 

-4 

200 

210 

1 

46438 

0 

8 

2225 

-215 

;  1 

469,40 

1  0 

-1 

-450 

240 

1 

,46942 

0 

+7 

-1300 

-70 

1, 

46944 

'  0.2  ; 

2 

275 

255 

1 

46945 

! 

0.2 

-7  , 

600 

195 

1  1/2 

46946 

-0.1 

+3 

1180 

-82 

1  1/2 

46948 

-0.1 

,  '  7 

-1200 

195 

2 

1 

46949 

-0.2 

11 

,-2000 

-125 

1  1/2 

46953 

! 

0.1 

»  1 

“2 

-1680 

-240 

2 

46956 

,  0.2 

1  20 

-25 

377 

2 

46958 

i 

-0.2 

8  , 

-833 

183 

I 
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8 .  CONCLUSIONS  AND  SUGGESTIONS 
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8.1  Summary  and  Conclusions 

(1)  Piezo  method  of  testing  gives  better  information  re- 
garding  pressures  inside  the  barrel  than  copper  crush¬ 
er  method  of  testing.  Piezo  gives  the  maximum  pressure 
at  the  gage  location.  It  also  gives  ignition  delay 
and  slope  which  are  indicators  of  propellant  character¬ 
istics  . 

(2)  Copper  crusher  deformation  can  be  considered  as  a 
weighted  integral  of  piezo  pressure-time  curve. 

(3)  Different  responses  from  piezo  and  copper  crusher 
are  correlated.  The  correlative  structure  is  given 
in  Section  2.5.  Correlation  analysis  is  found  to 
help  present  a  better  visual  picture  by  effecting 
data  reduction.  Copper  crusher  and  piezo  data  are 
found  to  exhibit  similar  pressure-velocity  relation¬ 
ships,  however,  copper  crusher  has  larger  variability. 
Piezo  data  indicate  a  time  trend  in  the  coating  pro¬ 
cess  which  crusher  data  fail  to  show.  This  time  trend 
could  be  significant  from  process  control  viewpoint. 

(4)  Acceptance  test  data  are  found  to  contain  time  trends, 
in  particular,  there  seems  to  exist  a  'point  of  shift* 
after  which  the  data  become  relatively  stationary. 

(5)  Method  of  cumulative  sum,  with  reference  constant 
equal  to  the  mean  of  the  series  of  data,  is  found  to 
give  the  best  visual  picture  of  the  point  of  shift. 


*-  asKw^wjs^^r^ 
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(6)  Data  after  the  point  of  shift  are  found  to  be  non¬ 
stationary  and  therefore,  nonnormal.  They  are  also 
autocorrelated . 

(7)  Analysis  cf  chamber  pressure  data  from  different 

manufacturers  indicates  the  mean  chamber  pressure 

to  be  quite  alike.  The  chamber  pressure  variance 

varies  considerably,  from  118  x  10  psi  for  Reming- 
4  T 

ton  to  S29  x  10  psi'  for  Winchester.  Mean  chamber 
pressure  of  Ball  ammunition  (46600  psi)  is  lower 
than  the  mean  chamber  pressure  for  Tracer  ammunition 
(49200  psi) .  The  within  lot  chamber  pressure  var¬ 
iance  for  Twin  City  is  29  x  10 5  psi2  and  for  Lake 

c  o 

City  is  17.6  x  10  psi  .  In  both  cases,  the  trend 
is  toward  a  reduction  in  variance,  indicating  con¬ 
tinued  improvement  in  production  and  testing  proes  ses, 

(8)  The  cutoff  date  and  the  corresponding  control  limits 
for  the  data  from  five  manufacturers  are  given  in 
Tables  4  and  5.  respectively.  On  the  average,  the  new 
control  .limits  are  found  to  be  lower  than  the  standard 
control  limit  of  52000  psi,  by  about  2000  psi. 

(S)  It  is  found  that  different  lot  selection  can  result 
if  ai.muniv.ion  selection  for  weapon  tests  is  based  on 
chamber  pressure  alone  or  port  p^ssure  alone.  Lot 
selection  is  therefore  based  on  the  bivariate  histo¬ 
gram  of  chamber  pressur.  and  port  p  essure. 
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(10)  The  selected  ammunition  lots  from  Lake  City  are 
given  in  Section  5.5  and  the  selected  lots  from 
Twin  City  are  given  in  Section  5.6, 

(11)  acceptance  testing,  standard  deviation  is 
determined  from  twenty  tests.  It  is  found  that 
twenty  tests  are  too  less  for  a  proper  estimate 
of  experimental  error. 

(12)  Propellant  lot  characteristics  are  found  to  be 
serially  correlated. 

(13)  Comparison  of  ammunition  data  from  B.A.A.P.  and 
loading  plants  indicate  no  change  in  propellant 
characteristics  during  the  elapsed  period  between 
the  tests. 

(.14)  The  model  building  approach  of  Appendix  I  appears 

to  be  the  correct  way  to  interpret  interior  ballis¬ 
tic  phenomena. 

8.2  Recommendations  for  Future  Work 

(1)  To  compare  copper  crusher  and  piezo  methods  as  a 

means  of  standard  testing,  following  analysis  should 
be  conducted. 

(a)  Cost  analysis  for  the  two  methods  should  be 
carried  out.  This  would  involve  the  deter¬ 
mination  of  proper  number  of  tests  to  be 
conducted  by  each  method.  Approach  of 
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Chapter  6  can  be  usjed  for  this  purpose. 

(b.)  An  analytical  expression  for  copper  crushes 

deformation  should  be  obtained  to  see  if  it 

> 

is  a  good  estimate  of  impact  energy. 

r 

(c)  From  the  viewpoint  of  gun  design,  efforts 

should  be  made  to  obtain  a  suitable  weighted 
integral  of  pressure-time  curve  as  a  measure 
of  pressure  inside  the  barrel.  Such  a  mea¬ 
sure  could  then  replace  the  present  measure 
of  crusher  deformation. 

(2)  Designed  experiments,  as  suggested  in  Section  5.7, 
should  be  carried  out  to  determine  the  ammunition 
factors  that  control  gun  performance.  These  factors 
can  be  used  to  refine  ammunition  selection  criteria. 

(3)  Analysis  in  Chapter  6  indicates  that  20  tests  are  not 
sufficient  for  proper  estimation  of  experimental  error. 
Further  analysis  should  be  conducted  to  determine  the 
proper  number  of  tests  necessary  for  acceptance  testing. 

(4)  Extensive  data  analysis  should  be  undertaken  to  deter¬ 
mine  if  identical  results  ire  obtained  at  P.A.A.P. 
and  Loading  plants.  For  this  purpose,  it  would  also 
be  necessary  to  evaluate  the  experimental  facilities 
at  different  loading  plants  and  B.A.A.P.  If  the 
results  are  identical,  then  there  is  a  possibility 

of  reduction  in  the  extent  of  acceptance  testing. 


(5)  Efforts  directed  at  building  a  model  for  interior 
ballistic  of  guns  are  likely  to  be  fruitful.  Sue! 
a  model  would  help  explain  the  observed  ballistic 
phenomena  and  suggest  likely  changes  for  improve¬ 


ment  in  the  weapon  system. 
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APPENDIX  I 


EMPIRICAL-MECHANISTIC  MODEL  FOR  INTERIOR  BALLISTICS  OF  GUNS 


Introduction. 


Over  the  past  several  decades,  considerable  effort  has 
been  directed  toward  the  theoretical  analysis  of  interior 
ballistics  of  guns.  The  complexity  of  the  phenomenon  has 
made  a  complete  theoretical  analysis  impractical.  Develop¬ 
ment  of  interior  ballistic  models  in  the  past  was  hampered 
by  lack  of  analytical  as  well  as  experimental  methods  and 
computer  facilities.  As  a  consequence,  several  simplifying 
assumptions  were  made.  In  view  of  the  current  emphasis  on 
high  velocity  weapons  and  the  considerable  improvements  in 
the  field  of  solid  propellants  {ball  powder,  improved  deter¬ 
rent  coatings,  etc.) ,  the  solutions  obtained  are  no  longer 
adequate.  With  today's  computer  facilities  and  the  recent 
advances  in  hydrodynamics  and  statistical  model  building 
techniques,  improved  interior  ballistic  models  appear  feasi¬ 
ble. 


Models  for  Physical  Systems 

Models  for  physical  systems  can  be  classified  into  three 
broad  categories  -  theoretical,  mechanistic  and  empirical. 
Theoretical  models  are  based  entirely  on  theoretical  consider¬ 
ations  and  presently  do  not  seem  feasible  for  interior  ballis¬ 
tics.  Mechanistic  models  assume  a  considerable  knowledge 
about  the  system  so  that  the  functional  form  of  the  model 
can  be  derived.  The  parameters  of  the  model  are  then 


estimated  from  the  data.  Empirical  models  are  not  based 
upon  a  physical  understanding  of  the  process.  Here  a 
functional  form  is  devised  and  parameters  estimated  to 
explain  the  observed  data.  For  the  interior  ballistics 
system  an  empirical-mechanistic  approach  seems  necessary, 
since  only  parted  information  is  available  about  the  system. 

Empirical-Mechanistic  Model 

To  build  a  mechanistic  model  there  is  a  need  to  know 
'how  the  system  works'.  Since  exact  information  about  the 
system  is  not  available,  true  functional  form  of  the  model 
cannot  be  determined.  However,  based  upon  partial  informa¬ 
tion  several  models  can  be  proposed  and  statistical  techni¬ 
ques  used  to  identify  the  most  probable  functional  form  of 
the  model.  From  experimental  results  the  parameters  of  the 
model  can  be  estimated .  The  model  is  then  diagnostically 
checked  for  adequacy  of  fit.  This  three  stage  iterative 
procedure  of  identification,  estimation  and  diagnostic  check¬ 
ing  will  lead  to  an  adequate  empirical-mechanistic  model . 

Once  the  model  is  available  it  can  be  used  to  predict 
the  process  behavior  under  different  conditions  or  for  opti¬ 
mization  of  the  process.  These  aims  can  also  be  achieved 
entirely  experimentally  using  empirical  models  and  response 
surface  techniques.  For  example,  optimum  burning  rate  can 
be  determined  by  making  powders  of  different  composition 
(by  varying  process  variables)  and  the  experimental  results 
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analyzed  using  response  surface  methods.  But  the  experi¬ 
mentation  involved  would  be  rather  expensive. 

Apart  from  optimization,  the  important  advantage  of 
this  approach  via  mechanistic  models  is  that  the  model  is 
useful  in  the  development  of  new  processes  since  meaningful 
extrapolation  is  possible.  Thus  the  ballistic  performance 
of  radically  different  experimental  weapons  or  powder  can 
be  simulated  using  these  models.  The  model  would,  therefore, 
be  of  considerable  use  in  developing  new  concepts  in  weaponary. 

Mathematical  Efodel  for  the  Gun 

As  a  first  step  a  simplified  model  is  assumed.  This 
is  based  upon  the  following  assumptions. 

A.  Assumptions 

(1)  Products  of  combustion  are  assumed  to  obey  ideal 
gas  law. 

(2)  Gases  of  combustion  are  assumed  to  have  negligible 
viscosity  and  thermal  conductivity.  Boundary  layers 
are  assumed  absent,  so  the  flov?  is  one  dimensional. 

(3)  Coulomb  friction  between  projectile  and  barrel  as 
well  as  the  resistance  due  to  compression  of  air 
ahead  of  the  projectile  is  neglected. 

(4)  Wal's  of  the  gun  are  assumed  perfectly  insulated, 
so  there  is  no  energy  loss  by  heat  transfer. 

Effect  of  chambrage  is  neglected. 


(5) 
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B.  Statement  of  the  Simplified  Problem 

In  short  the  problem  considered  is  as  follows:  There 
is  an  insulated  tube  sealed  at  one  end  and  open  at  the 
other  end.  Gases  and  powder  are  trapped  between  the  sealed 
end  and  the  pis ton ,  and  there  are  heat,  mass  and  volume 
sources  inside  the  volume  behind  the  piston  (Fig.  1).  The 
motion  of  the  piston  and  the  pressure  changes  constitute 
the  objects  of  the  solution  to  this  problem. 


Barrel  Projectile 


Figure  1 

C ’  Method  of  Solution 

The  Lagrangean  approach  .is  used  for  the  solution  of  this 
problem.  That  is  to  each  fluid  particle,  the  Lagrangean 
coordinate  x  is  assigned  to  be  its  distance  from  the  sealed 
end  at  time  time  t  -  0.  Then  at  any  instant  of  time,  t, 
its  location  is  given  by  the  function  R{x,  t) .  The  funda¬ 
mental  equations  governing  the  motion  of  gases  are  the 
equations  (1)  through  (5) . 
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Fluid  Equations: 

/ i \  3R  _ 

(1>  3t  *  U 


(2)  ir-  =  -V 

3t  vo  9x 


,,,  9E  _  9V  ,  * 

(3)-  at  "  aE  +  Q 


(“)  v  =  v  — 

°  3x 


(5)  E  =  f(p,V)  = 

Y-l 


Velocity  iquation 
Momentum  Equation 


Energy  Equation 
Continuity  Equation 


Equation  of  state. 


Here,  V,  u,  p,  E  are  the  specific  volume,  the  fluid  velocity, 
the  pressure  and  the  specific  energy  respectively.  Q  is  the 
rate  of  liberation  of  energy  and  y  is  the  ratio  of  specific 
heats,  VQ  is  the  reference  specific  volume.  These  equations 
can  be  directly  used  in  the  case  of  instantaneous  burning 
but  they  used  certain  modifications  for  finite  burning  rate. 

These  equations  of  fluid  motion  are  combined  with  the 
equation  of  motion  of  the  projectile,  given  by 


(6) 


£2  .  p  I 
,  ^ 

va  u  m 


where  m  is  the  mass  of  projectile  and  A  is  the  cross- 
sactional  area. 
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D.  Numerical  Integration  of  the  Set  of  Differential  Equations 

Since  an  analytical  solution  of  the  system  of  differen¬ 
tial  equations  is  impossible,  a  numerical  procedure  is  em¬ 
ployed.  The  scheme  is  taken  from  'Difference  methods  for 
initial  value  problems'  by  Richtmyer  and  Morton.  Certain 
modifications  have  been  made  to  take  care  of  the  finite 
burning  rate. 


Description  of  the  Experimental  Data  Used 

A  piezo  gauge  was  mounted  at  the  mid  chamber  position 
(point  A  in  Fig.  2(a))  and  the  pressure-time  history  was 
recorded  for  various  rounds.  A  typical  pressure-time  curve 
is  shown  in  Fig.  2(b).  Each  time  the  terminal  velocity  of 
the  projectile  was  measured  by  clocking  the  time  taken  by 
the  projectile  to  travel  between  two  magnetic  screens  placed 
at  a  known  distance  from  one  another. 


! - v  X 

i  A  Projectile 

-  -  V 

Chamber  Barrel 

.  2(a) 


Fig 


Fig.  2(b) 
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In  terms  of  the  experimental  data,  the' present  model 
building  problem  is  to  determine  an  appropriate  functional  ! 
form  of  Q  so  that  the  calculated  pressure-time  curve  and 
terminal  velocity  'match  as  closely  as  possible'  with  the 

I 

corresponding  experimentally  observed  responses,  i 

I 

Iterative  Model  Building  Procedure 

It  was  first  assumed  that  all  the  explosive  burns  before 
the  projectile  starts  to  move  (instantaneous  burning)  .  The'  ' 
resulting  pressure-time  history  is  shown  in  Fig.  3,  curve 
(i) .  Obviously  this  model  does  not  explain  the  experimental 
data.  As  a  matter  of  fact  the  curve  (i)  suggests  a  model 
for  Q  such  that  only  a  portion  of  the  explosive  burns  before 

i 

the  projectile  starts  to  move,  and  that  the  rest  of  the  powder 
continues  to  burn  while  the  projectile  is  being  accelerated. 
This  will  then  avoid  the  excessive  pressure  build  up  in  the 
chamber.  So  the  next  simple  model  was  to  assume  that  the 
powder  burns  at  some  constant  burning  rate  until  all  the 
powder  has  burned.  The  resulting  pressure  time  curve  is 
shown  in  Fig.  3,  curve  (ii) .  Although  it  is  a  substantial 
improvement  over  the  first  model,  it  still  does  not  explain 
the  experimental  curve  adequately.  The  main  differences  : 
being  that  the  model  exhibits  a  sharp  peak  wheras  the  experi¬ 
mental  curve  shows  a  smooth  variation  of  pressure  near  the 
peak  pressure  and  that  the  model  predicts  toe  high  a  pressure 
in  the  initial  portion  of  the  curve.  Requirement  of  lower 


I 
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pressures  in  the  beginning  suggests  the  use  of  slower  burn¬ 
ing  -ate  to  start  with  and  ah  increasing  burning  rate  as 
time  proceeds. 

A  study  of  the  burning  characteristics  of  the  explosives, 
in  the  mean  time,  revealed  that  the  burning  rate  should  be 
proportional  to  the  surrounding  pressure;  and  this  looked  to 

i 

be  an  ideal  functional  form  for  Q  because  initially  the  pres¬ 
sures  are  low  so  a  smaller  burning  rate  would  be  obtained  and 
as  the i pressure  starts  building  up,  the  burning  rate  would 
become  higher.  The  resulting  pressure-time  curve  is  shown 
in  Pig.  3,  curve  (iii)  .  It  is  apparent  by  observing  the 
similarities  in  this  curve  and  the  experimental  curve  that 

i 

by  suitably  estimating  the  proportionality  constant,  the 

experimental  curve  can  be  explained  except  for  the  portion 

i 

near  the  peak.  The  discrepancy  shows  that  Q  should  be  such 
that  after  a  portion  of  the  powder  has  burned,  the  burning 
rate  should  reduce.  This  was  then  taken  care  of  in  the  next 
model.  It  was  known  that  the  powder  is  used  in  the  form  of 
sperical  granulations.  So  as  the  burning  proceeds,  the  sur¬ 
face  ar^a  of  the  powder  should  go  on  reducing.  It  was  also 
known  that  a  deterrent  coating  is  applied  on  the  powder  and 

the  effect  of  this  was  assumed  to  nullify  the  effect  of 

' 

change  surface  area.  So  the  new  model  for  Q  is 

,  i 

Q  =  HF  x  RMBC  x  P,  until  (1-a)  x  100%  of  powder  has  burned 
=  HP  x  RMBC’  x  P  x  surface  area,  for  the  remaining  portion 


of  mass. 
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RMBC  and  RMBC*  are  such  that  the  two  equations  yield 
ti\e  sajne  q  when  transition  takes  place.  HF  and  P  are  res¬ 
pectively  uhe  calorific  value  and  the  pressure.  The  result¬ 
ing  pressure  time  curve  is  given  in  Pig.  3,  curve  (iv) .  It 
is  apparent  that  this  curve  explains  the  basic  nature  of  the 
experimental  curve.  Now  it  was  thought  worthwhile  to  esti¬ 
mate  the  parameters  a,  RMBC,  and  the  Primer  energy  (taken 
as  the  small  amount  of  energy  instantaneously  available  at 
the  start  of  combustion)  to  get  the  best  possible  fit.  In 

Fig.  4  the  five  available  experimental  curves  are  sketched. 

2 

On  the  assumption  that  the  errors  are  NID(Of  Q  ),  least 
squares  criterion  is  used  and  parameters  estimated  using 
UWHAUS.  The  fitted  curve  is  also  sketched  on  Fig.  4. 

(Note:  the  data  was  discretized  by  taking  41  points  from 
each  curve) . 

Using  these  estimated  values  of  parameters,  pressure  and 
velocity  histories  along  the  length  of  the  barrel  were  deter¬ 
mined.  Figure  5  shows  the  pressure  variation  along  the 
length  of  the  barrel.  Figure  6  shows  the  breach  pressure 
and  projectile  base  pressure  variation  as  a  function  of  the 
distance  travelled  by  the  projectile.  Breach  pressure  is 
always  found  to  be  higher  than  the  projectile  base  pressure. 
The  figure  also  shov7S  the  variation  of  projectile  velocity 
along  the  length  of  the  barrel.  The  predicted  muzzle  velocity 
is  found  to  be  much  lower  than  the  actual  muzzle  velocity. 

The  variation  of  breach  pressure  and  projectile  base  pressure 
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as  a  function  of  time  is  indicated  in  Fig.  7.  This  breach 
pressure  curve  is  the  one  that  is  being  compared  with  the 
piezo  pressure-time  curve.  Finally,  the  variation  of  pro¬ 
jectile  velocity  and  distance  with  time  is  shown  in  Fig,  8, 

Results 

In  appropriate  units,  the  parameter  values  and  their 
confidence  intervals  are 


a  =  0.3491 


PRIMER  =  .1504  x  10 


>10. 


-11 


RMBC  =  2.385 


a:  (.2917  to  .4064),  PRIMER  x  10  :  (-.84  to  .87) 

RMBC  :  (2.378  to  2.392) 


Correlation  matrix  a 

1.0 


PRIMER 

-.38 

1.0 


RMBC 

-.19 

.06 

1.0 


Analysis  of  Variance 

#  of  obs,  =  41  x  5  =  202,  #  of  parameters  =  3 

Residual  sum  sq.  =<  .27845  x  10*^,  d,  f.  =  202 

16 

Pure  error  sum  sq.  =  .02649  x  10  ,  d.  f.  =  4  x  41  =  164 
Lack  of  fit  sum  sq,  =  0.24856  x  10^,  d.  f.  =  38 
Therefore, 

F  statistic  -  ^mjLiSlZ3S_ 

.02649  x  10l6/164 


=  40.6  (38,164  degrees 
of  freedom) 
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Looking  into  F  tables,  for  10%  significance  level  the 
corresponding  F  value  is  approximately  1.30.  The  model  is 
found  to  be  inadequate. 

Analysis  of  Residuals,  Diagnostic  Checking 

A  look  at  the  residuals  tells  that  the  model  is  inade¬ 
quate,  as  also  seen  from  the  r  test.  Additionally,  the  trend 
in  residuals  shows  where  the  inadequacy  lies  and  how  to  im¬ 
prove  upon  it. 

(a)  The  fitted  values  are  consistently  lower  in  the 
initial  part  of  the  curve.  This  may  mean 

(1)  higher  initial  burning  rate 

(2)  increased  engraving  force 

(3)  increased  primer  energy 

(b)  By  bodily  translating  the  experimental  curves  to 
the  right,  a  better  fit  for  the  initial  part  of 
the  curve  (i.e.  the  portion  during  which  the  pres¬ 
sure  is  increasing)  is  obtained.  This  has  a  rele¬ 
vance  to  the  experimental  conditions  because  of 
the  possibility  of  some  time  lag. 

(c)  The  fitted  values  have  a  large  curvature  near  the 
peak  pressure  as  compared  with  that  of  the  experi¬ 
mental  data.  A  smaller  burning  rate  near  this 
point  would  remedy  this  situation. 

(d)  The  effect  of  the  assumption  of  lack  of  heat  trans¬ 
fer  and  the  energy  equation  used  should  be  examined 
further . 


(e)  The  above  modifications  are  expected  to  also 
take  care  of  the  one  sided  residuals  in  the  • 
position  of  the  curve  where  the  pressure  is 
dropping . 

Note  on  Estimation  Procedure 

It  is  more  realistic  to  take  the  errors  to  be  correlated 
when  data  is  obtained  from  one  experimental  curve.  However, 
from  one  experimental  curve  co  the  next  experimental  curve, 
the  errors  can  be  assimed  uncorrelated.  Then  one  can  think 
of  an  estimation  criterion  similar  to  the  one  described  by 
Box  and  Draper  in  their  paper  "The  Bayesian  Estimation  of 
common  parameters  from  several  responses,"  which  will  be  more 
appropriate  for  the  present  problem. 

The  velocity  has  been  totally  ignored  in  the  estimation 
procedure  followed  so  far.  This  is  because  when  the  estima¬ 
tion  is  done  with  the  pressure  response  alone,  the  expected 
velocity  corresponding  to  the  best  fit  is  about  15%  smaller 
than  the  observed  velocity.  Hence  if  the  estimation  is  carried 
out  with  only  the  response  velocity,  a  disjoint  confidence 
region  for  the  parameters  would  result.  So  it  does  not  seem 
necessary  to  carry  out  a  multiresponse  analysis  with  the 
present  model. 

The  confidence  interval  of  the  parameter  PRIMER  does 
include  the  point  zero.  However,  because  of  some  physical 
conditions,  it  was  chosen  to  retain  it  in  the  estimation 
problem. 
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*  .  4  * 

Uses  of  the  Empirical-Mechanistic  Model 

i  • 

The  fitted  empirical-mechanistic  model  can  be  usv.u  in 
several  ways.  Some  of  the  uses  are  based  upon  the  fact  that 

i 

the  model  explains  observed  data.'  Others  stem  from  the  me¬ 
chanistic  aspects  of  the  model  that  permit  meaningful  extra- 

l 

polation.  ’ 

(1)  The  model  can  be  used  to  simulate  the  interior  ballistic 
performance  of  a  weapon  for  any  given,  set  of  initial 
conditions.  This  is  especially  useful  in  situations 
where  actual  experimentation  is  difficult,  as  in  eva- 

l 

luating  the  effect  of  holding  force. 

(2)  The  model  is  useful  for  process  optimization.  For  exam¬ 
ple,  burning  characteristics , can  be  optimized  to  yield 
the  desired  muzzle  velocity  with  small  internal  pressures. 
An  anticipated  problem  in  this  regard  is  the .development' 
of  a  suitable  objective  functibn  for  internal  pressures 
to  be  minimized. 

1  '  i 

(3)  Once  process  optimum  is  defined,  it  can  be  attained  by 

I 

suitable  changes  in  process  variables.  Because  of  the 
disturbances  entering  the  system,  the  process  would 

.  i 

continually  deviate  from  the  optimum  necessitating  an 

I 

efficient  testing  procedure.  Such  a  testing  procedure 
should  detect  changes  quickly  and  also  point  out  ^ssignr 
able  causes.  i 

The  empirical-mechanistic  model  can  be  used  to  establish 
such  a  testing  procedure.  The  effect  of,  different 


I 


aS2 

.  i'  ■  . 

i 

process  variables  on  (say)  piezo  pressure-time  curve 
1  can  be  .simulated  using  this  model.  The  observed 

i  ^ 

,  preSsure-time  curve  would  then  indicate  the  process 
variables  .that  need  to  be,  adjusted.  It  is  possible 
that  several  process  variables  might  influence  the 

i 

piezo  curve’  in  the  same  fashion  thereby  making  detect¬ 
ion  difficult.  This  difficulty  might  be  circumv anted 
by  considering  multiple  responses  that  distinguish 
between  the  effects  of  different  process  varieties. 

(4)  The  empirical -mechanistic  model  can  be  used  to  make 
small  scale  replicas  of  large  weapons,  thereby  reducing 

i 

cost  of  weapon  testing. 

(5)  The  model  can  be  used  to  evaluate  the  performance  of 

I 

radically  different  experimental  weapons  or  powders. 

1  For  example,  the  effect  of  encapsulation  or  of  propel¬ 
lants  with  lower  molecular  weights  or  of  multiple  charges  . 
can  be  studied  with  this  model.  The  model  can  be  used 
to  determine  process  changes  that  would  lead  to  increased 
velocity.  Such  inferences  are  possible  only  because 

I 

the  model  is  based  upon  mechanistic  considerations. 
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COPPER  CRUSHER  AND  PIEZO  DATA  FROM  FRANKFORD  AND  B.A.  A.  P. 
FRANKFORD  ARSENAL  TEST  DATA 

*.>  rmmir^^nwaM— — mi  iw— — u— 


Special  Ammunition  Tests 


Cartridge  -  5.56mm  (Ball) 

o 

Temp.  -  70  F 


Pressure 


Round 

# 

Peak  Port 
Pressure 
(xlOOpsi) 

Peak  Chamber 
Pressure 
(x  lOOps i) 

Velocity 

(fps) 

Slope 

Time 

Integral 

(psi-sec) 

Action 

Time 

(MS) 

1 

X 

130 

465 

3161 

3.376 

27.93 

1.16 

Lot 

LC12604 

2 

130 

490 

3207 

4.165 

27.40 

1.237 

(High) 

3 

125 

465 

3158 

3.376 

28.79 

1.191 

4 

130 

500 

3192 

3.606 

29.43 

1.148 

5 

125 

485 

3181 

3.487 

29.40 

1.163 

6 

125 

495 

3203 

3.732 

28.98 

1.261 

7 

130 

485 

3177 

3,376 

29.18 

1.235 

8 

130 

475 

3180 

3.487 

29.04 

1.29 

9 

120 

485 

3176 

3.376 

29.20 

1.164 

10 

130 

495 

3193 

3.606 

29.23 

1.265 

Lot 

LC12594 

(Low) 

1 

120 

485 

3172 

3.172 

29.92 

1.238 

2 

125 

505 

3193 

3.271 

28.23 

1.225 

3 

120 

590 

3237 

3.487 

28.70 

1.184 

4 

120 

495 

3172 

3.271 

28.98 

1.227 

5 

120 

520 

3221 

3.078 

29.07 

1.184 

6 

125 

510 

3222 

3.867 

28.57 

1.155 

7 

125 

510 

3185 

3.078 

29.26 

1.152 

8 

125 

515 

3237 

3,867 

28.96 

1.165 

9 

125 

515 

3223 

3.487 

28. 90 

1.316 

10 

125 

515 

3225 

3.376 

28.68 

1.244 

! 

r 
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APPENDIX  II  (B) 

BADGER  DATA  (COPPER  CRUSHER) 

COMPOSITES 
COMPONENT  -  TWIN 
CALIBER  -  5.56  MM  (BALL) 

TEMPERATURE  =  70°F 


Coating 

Date 

Composites 

Chamber 

Pressure 

(xlOOpsi) 

Velocity 

from 

Velocity 

Barrel 

(fps) 

Velocity 

from 

Pressure 

Barrel 

(fps) 

Charge 

Weight 

(Grains) 

7-  3-70 

244 

429 

3251 

3204 

27.9 

7-  7-70 

241 

517 

3245 

3226 

27.6 

7-  9-70 

243 

469 

3258 

3230 

27.7 

7-  9-70 

246 

499 

3250 

3240 

25.8 

7-  9-70 

245 

450 

3243 

3234 

27.6 

7-10-70 

242 

424 

3240 

3192 

26.1 

7-14-70 

248 

469 

3255 

3235 

27.5 

7-14-70 

250 

524 

3246 

3234 

27.3 

7-14-70 

252 

479 

3254 

3234 

27.7 

7-17-70 

251 

472 

3252 

3273 

26.5 

7-17-70 

253 

446 

3245 

3226 

27.5 

7-21-70 

249 

447 

3250 

3250 

27.4 

7-21-70 

258 

416 

3254 

3158 

26.4 

7-21-70 

259 

459 

3250 

3192 

28.0 

7-21-70 

254 

432 

3253 

3182 

27.5 

7-21-70 

255 

479 

3246 

3189 

27.6 

7-21-70 

256 

420 

3239 

3179 

27.9 

7-22-70 

260 

438 

3255 

3215 

28.1 

7-22-70 

261 

455 

3241 

3216 

28.3 

7-27-70 

263 

453 

3245 

3200 

28.3 

7-27-70 

265 

435 

3247 

3187 

27.7 

7-30-70 

264 

483 

3253 

3213 

26.5 

7-31-70 

267 

482 

3249 

3221 

28.0 

S-  4-70 

266 

479 

3254 

3272 

28.0 

8-  4-70 

268 

489 

3254 

3216 

26.5 

185 
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Coating 

Date  Coapcaites 

Chamber 

Pressure 

(xlOOpsi) 

Velocity 

from 

Velocity 

Barrel 

(fps) 

Velocity 

from 

Pressure 

Barrel 

(fps) 

Charge 

Weight 

(Grains) 

8-  4-70 

269 

449 

3239 

3201 

26.5 

8-  4-70 

270 

491 

3250 

3213 

28.0 

8-  4-70 

271 

461 

3240 

3230 

28.3 

8-  7-70 

272 

461 

3254 

3210 

28.1 

8-  7-70 

273 

447 

3252 

3199 

27.9 

8-  7-70 

274 

439 

3242 

3198 

26.5 

8-  7-70 

275 

419 

3255 

3192 

27.9 

8-  7-70 

277 

485 

3254 

3235 

27.8 

3-10-70 

276 

460 

3251 

3202 

28.0 

8-10-70 

278 

499 

3246 

3200 

28.2 

8-13-70 

279 

440 

3247 

3189 

27.7 

8-13-70 

280 

467 

3247 

3182 

26.2 

8-13-70 

281 

488 

3258 

3202 

28.4 

8-17-70 

283 

420 

3244 

3166 

28.0 

8-17-70 

284 

478 

3242 

3172 

27.4 

8-17-70 

285 

425 

3250 

3170 

28.4 

8-21-70 

286 

439 

3241 

3174 

26.3 

8-21-70 

288 

446 

3246 

3223 

27.9 

8-24-70 

289 

455 

3250 

3206 

25.9 

8-24-70 

290 

428 

3244 

3214 

27.8 

8-24-70 

291 

467 

3249 

3227 

27.3 

8-24-70 

292 

447 

3254 

3214 

27.  S 

8-26-70 

293 

448 

3242 

3222 

27.5 

8-26-70 

294 

499 

3246 

3231 

26.0 

8-31-70 

296 

441 

3252 

3233 

28,0 

8-31-70 

298 

449 

3249 

3216 

27.8 

9-  3-70 

297 

467 

3241 

3222 

27.6 

9-  3-70 

300 

487 

3246 

3235 

27.5 

9-  4-70 

301 

474 

3249 

3730 

26.2 

9-  4-70 

302 

444 

3252 

3328 

27.8 

9-  4-70 

303 

459 

3243 

3240 

37.9 

186 


Coating 

Date  Casapoeitcs. 


9-  4-70 

304  „ 

9  -  4-70 

308 

9-  9-70 

309 

9-  9-70 

314 

9-11-70 

306 

9-14-70 

307 

9-14-70 

310 

9-14-70 

313 

9-23-70 

315 

9-23-70 

316 

9-23-70 

324 

9-24-70 

320 

9-24-70 

326 

9-25-70 

317 

9-25-70 

313 

3-25-70 

335 

Velocity 
>  i-from . 


Chamber 

Pressure 

(xlOOpsi) 

Velocit 

Barrel 

(fps) 

n  477 

3242. 

■'  432 

3245 

477 

3240 

444 

3245 

473 

3256 

481 

3240 

483 

3253 

427 

3249 

457 

3240 

432 

3256 

443 

3245 

453 

3251 

451 

3250 

538 

3252 

503 

3252 

462 

3238 

Velocity 

from 

Velocity  Charge 

-  Barrel  Weight 

(fps)  (Grains) 


3242 

27.8 

3234 

27.8 

3233 

27.1 

3212 

27.5 

3240  - 

26.0 

3196 

27.3 

3198 

27.6 

3179 

27.5 

3199 

28.2 

3204 

25.9 

3213 

27.6 

3241 

28.5 

3234 

27.9 

3246 

28.1 

3230 

27.6 

32224 

26.1 
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APPENDIX  II  (C) 

BADGER  DATA  (Pieso  Transducer) 

.  t  ;  .  -  V  .  .  ‘  ” 

,  i  Composites 

Component  -  Twin 
Caliber  -  5.56rtgn  (Ball) 

Charge  Weight  Sames  As  In  Copper  Crush 

(x  100  psi)  125° 


Coating 

date 

Composites 

Peak  Chamber 
Pressure 

Velocity^ 
(fps) 70° 

Velocity, 
(fps) 125° 

Peak  Chamber 
Pressure 

7-  3-70 

244 

504 

501 

7-  7-70 

241 

543 

3256 

3295 

550 

7-  S-70 

243 

504 

3243 

3272 

526 

7-  9-70 

246 

518 

3225 

3269 

542 

7-  9-70 

245 

501 

3225 

3245 

500 

7-10-70 

-  242 

498 

3209 

3267 

529 

7-14-70 

248 

521 

3248 

3282 

527 

7-14-70 

250 

535 

3253 

3275 

537 

7-14-70 

252 

527 

3261 

3295 

545 

7-17-70 

251 

509 

3229 

3285 

541 

7-17-70 

253 

505 

3213 

3256 

528 

7-21-70 

249 

507 

3246 

3284 

532 

7-21-70 

258 

503 

3223 

3283 

542 

7-21-70 

259 

506 

3226 

3272 

520 

7-21-70 

254 

515 

3238 

3271 

530 

7-21-70 

255 

524 

3252 

3283 

544 

7-21-70 

256 

526 

3265 

3269 

526 

7-22-70 

260 

520 

3239 

3281 

535 

7-22-70 

261 

505 

3229 

3288 

531 

7-27-70 

263 

510 

3227 

3284 

521 
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C^&tlrsgy, 

: V 

Ck  100  psi) 70° 

.  Peak  Chamber 
Pressure 

•,»  *  « r.-  •  „T  ‘  C  1 

Velocity, .. 
<fpa)7Q° 

**  t 

iVelocitv 

(fpsU25° 

(x  100  psi) 125° 
Peak  Chamber 
Pressure 

7-27- 70 

265 

498 

3206 

•>  r 

3274 

525 

7-30-70 

264 

494 

V  i  .  "j 

3225 

3260 

513 

7  31-70 

267 

516 

3216 

3274 

546 

8-  4-70 

256 

519 

3252 

3289 

i 

525 

8-  4-70 

268 

494 

3194 

3263 

531 

8—  4-70 

265 

505 

3204 

3303 

567 

8-  4-70 

270 

501 

3179 

3244 

529 

8-  4-70 

271 

488 

3206 

3267 

538 

8-  7-70 

272  • 

469 

3149 

3264 

533 

8-  7-70 

273 

453 

3182 

3245 

521 

8-  7-70 

274 

501 

3198 

3270 

539 

8-  7-70 

275 

492 

3198 

3259 

538 

8-  7-70 

277 

520 

3237 

3259 

538 

8-10-70 

276 

498 

3214 

3255 

521 

8-10-70 

278 

508 

3212 

3252 

528 

8-13-70 

279 

495 

3198 

3241 

515 

8-13-70 

280 

491 

3185 

3241 

521 

8-13-70 

281 

497 

3212 

3242 

513 

8-17-70 

283 

488 

3188 

3258 

522 

8-17-70 

284 

481 

3175 

3217 

504 

8-17-70 

285 

492 

3206 

3251 

510 

8-21-70 

286 

499 

3170 

3271 

590 

8-21-70 

288 

543 

3211 

3257 

571 

8-24-70 

290 

510 

3176 

3217 

551 

8-24-70 

291 

513 

3189 

3240 

558 

8-24-70 

292 

502 

3188 

3245 

566 

881 
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Ztl (leq  00 

iGcatifrig : 
Jtettae  ■' 

Jt  X' 

&'j‘i  ,  -  J  L'J‘ 

Composites 

<X  100  psi)70° 

■ '  i  o  /Peak  Chamber 
Pressure 

v  :■*  . 

velocity, 
(fps) 70  0 

Velocity, 
(fps) 125° 

(x  100  psi) 125 
Peak  Chamber 
Pressure 

8-24-70 

292 

515 

3207 

3225 

535. 

8-26-70 

293  • 

501 

3195 

3244 

543 

8-26-70 

294 

523 

3186 

3211 

549 

8-31-70 

296 

531 

3232 

3215 

543 

8-31-70 

298 

533 

3201 

3228 

551 

9-  3-70 

297 

538 

3243 

3252 

560 

9-  3-70 

300 

534 

3226 

3275 

578 

9-  4-70 

301 

551 

3242 

3286 

588 

9-  4-70 

302 

558 

3280 

3311 

569 

9-  4-70 

303 

557 

3266 

3300 

578 

9-  4-/0 

304 

548 

3262 

3288 

568 

9-  4-70 

308 

534 

3238 

3306 

582 

9-  9-70 

309 

529 

3215 

3281 

584 

9-  9-70 

314 

534 

3249 

3283 

576 

9-11-70 

306 

577 

3272 

3308 

599 

9-14-70 

307 

547 

3241 

3331 

600 

9-14-70 

310 

560 

3243 

3299 

584 

9-14-70 

313 

519 

3218 

3314 

566 

9-23-70 

315 

553 

3225 

3272 

570 

9-23-70 

316 

548 

3223 

3273 

597 

9-23-70 

324 

524 

3239 

3270 

565 

9-24-70 

320 

511 

3249 

3298 

577 

9-24-70 

326 

480 

3209 

3257 

560 

9-25-70 

317 

534 

3244 

3315 

595 

9-25-70 

319 

536 

3229 

3297 

585 

9-25-70 

335 

490 

2214 

3293 

578 
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APPENDIX  XX  (D)  . 

BADGER  DATA  (Piezo  Transducer) 
COMPOSITE 
COMPONENT  -  TWIN 
CALIBER  -5.56  (BALL) 

’  TEMPERATURE  -  70°F 


Composites 

Slope 

Ignition 
Delay  (ms) 

241 

6.08 

0.18 

242 

4.0 

0.30 

243 

5.31 

0.20 

244 

4.66 

0.28 

245 

4.08 

0.26 

246 

4.79 

0.27 

248 

4.24 

0.27 

249 

4.7 

0.22 

250 

5.78 

0.24 

251 

3.88 

0.32 

252 

5.08 

0.23 

253 

4.7 

0.24 

254 

4.92 

0.25 

255 

5.32 

0.23 

256 

4.74 

0.26 

258 

3.9 

0.29 

c.fr 1 


APPENDK  n  (E) 

BADGER  DATA  (COPPER, CRUSHER) 


191 


Ban<3  Blends 
For  Twin  City 


Caliber  -  5.56mm (Ball) 
Temp.  -  70°F 


Sample 

No. 

Coating 

Date 

Charge 

Weight 

(grains) 

Velocity 
(Vel.  Barrel) 
_Jfps) 

Chamber 

Pressure 

(psi) 

Velocity 

(Pressuj 

(fps) 

HB-644 

7-  2 

27.1 

3249 

459 

3232 

HB-657 

7-16 

27.1 

3249 

452 

3205 

HB-662 

7-24 

27.0 

3244 

460 

3187 

HB-663 

7-27 

27.5 

3241 

449 

3183 

HB-666 

7-28 

27.6 

3239 

480 

3211 

HB-674 

8-  4 

27.4 

3246 

447 

3199 

HB-678 

8-  6 

27.3 

3249 

460 

3208 

HB-680 

8-17 

28.1 

3242 

462 

3184 

HB-692 

8-25 

27.3 

3248 

478 

3239 

HB-697 

8-31 

27.4 

3251 

457 

3224 

HB-717 

9-18 

27.2 

3244 

460 

3219 

HB-719 

9-23 

27.6 

3248 

488 

3253 

HB-720 

9-23 

28.0 

3254 

465 

3230 

APPENDIX  n.(F). 

BADGER  CATA’ICOPPER  CRUSHER) 


i„j  , 
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Hand  Blende  *  c. 
For  Lake  City 


.  ,  *■»  « 

Caliber  - 

5. 56mm (Ball) 

■* 

A  ?>'  ^ 

Temp.  -  70 

°F 

Sample 

No. 

{  * 

Coating 

Date 

Charge 

Weight 

(grains). 

Velocity 
(Vel.  Barrel) 
(fps) 

Chamber 

Pressure 

(psi) 

Velocity 

(Pressure  Barrel) 
(fps) 

'  r 

HB-646 

7-  7 

28.0 

3240  ' 

440 

3239 

HB-647 

7-  7 

27.9 

3251 

472 

3232 

HE-648 

7-  8 

26.9 

3247 

456 

3214 

HB-652 

7-13 

27.5 

3242 

466 

3245 

HB-653 

7-13 

27.3 

3251 

476 

3250 

HB-661 

7-23 

27.6 

3240 

477 

3254 

HB-667 

7-30 

27.8 

3248 

464 

3209 

HB-679 

8-10 

28.2 

3248 

454 

3222 

HB-688 

8-21 

27 .4 

3249 

461 

3250 

HB-695 

8-27 

21  A 

3242 

443 

3226 

HB-701 

9-  3 

27.9 

3247 

477 

3232 

HB-713 

9-16 

27.1 

3238 

464 

3253 

HB-715 

9-17 

27.0 

3252 

469 

3241 

HB-719 

9-23 

27.6 

3248 

488 

3254 

HB-726 

9-28 

28.0 

3245 

455 

3256 

HB-730 

10-2 

27.4 

3245 

462 

3265 

; :  .-fi  w-*  '  -  ;•.»  .  Caliber  -  5.56  (Ball) 

~  .  Teiap.  -  70°P 

Charge  Weight  Same 

+  as  for  Copper  Crush 

!  1  3 i  j *  \  ?I  . 

-  Peak  ..Chamber  Pressurs-time 

Sample  Coating  Pressure  Velocity  Integral  Ignition 


No, 

Date 

(psi) 

_  (fps)  . 

(psi-sec) 

Slope 

(ms) 

KE-S44 

7-  2 

547 

3250 

27-773 

2.111 

0.23 

HB-657 

7-16 

540 

3232 

25-420 

2.179 

0.24 

HB-662 

7-24 

514 

3215 

HB-663 

7-27 

481 

3168 

X 

X 

X 

HB-666 

7-28 

514 

3224 

X 

X 

X 

HB-674 

8-  4 

483 

3180 

26-470 

1.625 

0.29 

HB-678 

8-  6 

489 

3172 

26-456 

1.516 

0.30 

HB-680 

8-17 

491 

3200 

26-940 

1.799 

0.25 

H B-692 

8-25 

525 

3220 

26-039 

2.050 

0.30 

HB-697 

8-31 

530 

3216 

26-287 

1.927 

0.30 

HB-717 

9-18 

511 

3237 

25-420 

1.958 

0.29 

HB-? 19 

9-23 

537 

3233 

26-955 

2.267 

0.28 

HB-720 

9-23 

555 

3268 

28-697 

2.117 

0.25 

APPENDIX  11(H) 
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BADGER  DATA  (PIE20  TRANSDUCER) 

Hand  Blends 
For  Lake  City 

V-  J  ,  *  ’  Caliber  ^  ;5.>56nun(iaii) 

Temp.  -  70°F 

t*  *  ’  * 

*  \  ' 

Charge -Weight  Same 
as  for  Copper  Crush 


Sample 

No. 

Coating 
Date  ' 

Peak  Chamber 

Pressure 

(psi). 

Velocity 

(fps) 

Pressure-time 

Integral 

(psi-sec) 

Slope 

Ignition 

Delay 

(ms) 

HB-646 

7-  7 

532 

3224 

28-167 

1.67 

0.19 

HB-647 

7-  7 

536 

3225 

27-877 

1.646 

0.19 

HB-648 

7-  8 

532 

3210 

29-307 

1.884 

0.19 

HB-652 

7-13 

542 

3243 

28-813 

1.828 

0.20 

HB-653 

7-13 

536 

3233 

28-467 

1.854 

0.18 

HB-661 

7-23 

509 

3245 

X 

X 

0.19 

HB-667 

7-30 

519 

3252 

X 

X 

0.16 

HB-679 

8-10 

520 

3264 

28-806 

1.754 

0.19 

HB-688 

8-21 

554 

3248 

28-827 

2.028 

0.21 

HB-695 

8-27 

522 

3216 

28-136 

1.896 

0.23 

HB-701 

9-  3 

532 

3256 

X 

X 

0.21 

HB-713 

9-16 

536 

3265 

27-819 

1.961 

0.20 

HB-715 

9-17 

539 

3265 

27-627 

2.002 

0.20 

HB-719 

9-23 

538 

3255 

28-67 

2.158 

0.25 

HB-726 

9-28 

465 

3121 

X 

X 

X 

HB-736 

10-2 

447 

3104 

X 

X 

X 

&p{ 

\Ht  l*  .<iQi-;  lifj> 

■ 

,  -  '  ,V> •','!£•  «>  '.  i’  \t~!  xIm  A’'.jV  .  ’ 

X95 

Appendix  in 

i 

Data  for  Sequential  Variance4  Estimation 

\*  */  ^  * 

■ 

Chamber,  .fereigsure  Data  froraBAAt*  Piezo  Data  f  rom  7~7, 

i  «WTW»^tea£)  «  Frankford  Arsenal (psi) 

m***""  ■  ■  "■'">'■1  . . ffin.i  !■■—  i  —  . . . .  —  . .  ■  "1,1.1  ■■'■■! 


Serial 

Number 

Ltffc'O** 

46881 

Lot . 
46893 

Lot 

46626  1 

LC12604 

LC12594 

1  "  >r 

4  (>400 

45000 

43700 

46500 

48500 

2 

445.00 

44300 

43200 

:  49000' 

i 50500  , 

3 

43600 

42900 

44900 

46500 

5<[000 

4 

45200 

44900 

44700 

50000 

•  49500. 

5 

45600 

44600 

43600 

■  48500 

52000 

6 

46600 

44700 

43800 

49500 

51000 

i  ' 

7 

45800 

45300 

45700  i 

48500 1 

51000 

8,. 

46400 

44500 

45100 

47500 

51500 

9 

44700 

45200 

44300 

48500 

51500 

10 

48100 

44700 

45900  : 

49500 

51500 

11“ 

46300 

42700 

45500 

l 

12  , 

44000 

44400 

43800 

. 

1 

13 

46100 

44200 

42800 

!  1 

14 

47100 

45700 

42800 

1 

\ 

15 

45400 

44900 

46700 

16 

45200 

42700 

45700 

t 

17 

47500 

43200 

44100 

* 

18 

43800 

43300 

45100 

19 

43800 

43400 

43800' 

20 

47000 

46100 

42800 

T 
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PROPELLANT  1  CHARGE  WEIGHT  AND  CHAMBER  I 


, —  ’  — PRESSURE-  DATA'  FROM -BAAP 

■  -  ■  - 

■*  fu- 

Date, 
Fired  -  - 

,  A  " 

*1 

.'■7  '<  ' 

Lot  No. 

J  '  'IT 

*  -  *  i 

-Component^ 

t  Y>  Charge 
j,*.  Weight , 
■(Grains).: 

Corrected  Mean 
Chamber  Pressure 
x  100  psi  . 

12-31-68 

1  j 

457  47 

1  LC 

26.9 

479 

1-10-69 

■'  45748 

TW 

26.fc 

475 

1-14-69 

45749 

?  ’*  i 

TW 

27.2 

490 

1-16-69 

4575,0 

LC 

27.3 

492 

! 1-16-69 

45750 

1  FA 

27.0 

495 

1-22-69 

45751 

TW 

27.2 

467 

1-27-69 

45,752 

TW 

‘  27.3 

479 

2-  7-69 

45753 

TW  ; 

27.2 

476 

1  2-.11-69 

46^94 

T¥f 

,  27.3  ' 

471 

2-13-69. 

i  46195  ’ 

LC 

27.2 

489 

2-17-69 

46196 

TW 

27.6 

482 

2-18-69 

•  46196 

LC 

27.9  « 

480 

2-21-69 

46198 

TW 

27.4 

467 

2t22-69 

46199 

LC 

27.6 

473 

2-24-69 

,  46200 

i  LC 

27.? 

482 

2-26-69 

46201 

TW 

27.3 

477 

2—28—69 , 

i 

46202 

TW 

27,3 

472 

3-  4-69 

46203 

LC 

27.1 

458 

3-  5-69 

46204 

TW  , 

27.2 

456 

3-  6-69 

46205 

TW 

27.0 

445 

3-12-69 

46206 

!lc 

27.0 

488 

3-17-6? 

46,207 

TW 

26.9 

445 

3-18-69 

462Q8 

TW 

26.9 

457 

3-20-69 

46209 

LC 

27.2 

455 

3-21-69 

<  46210 

J 

TW 

26.7  , 

449 

3-24-69 

46211 

TW 

27.1 

450 

3-25-69 

46212 

LC 

27.3 

469 

,3-27-69 

46213 

TW 

26.6 

462 

- - SMi  ta- jgj  - — 

- - ~. ':_: — charge-  -  -  Gorrv^tecbMean - 

Datfe^W  s  >  ’’Weight  Chamber:  Pressure 

Fitsasiibt^Niyi^CecaponehiB^  ‘(Grains)  x  100  psi  •-‘,pu 

.  Liaj  a.j.f  r  ■  . iiaa^.5lA.«<w  ,;. 


4-  7-69 

46214 

LC ' 

26.8 

471 

4-  7-69 

46215 

TW 

26.7 

455 

4-  9-69 

4625.6 

TW 

26.5 

482  ' 

4-10-69 

4621*7 

TW 

26.8 

481 

4-11-69 

46^11 

--P.S 

LC  " 

27.1 

480 

4-16-69 

46219 

TW 

26.8 

470 

4-17-69 

46221} 

LC 

27.3 

448' 

4-22-69 

46221 

TW 

27.0 

478 

4-24-69 

46222 

TW 

27.3 

463 

4-25-69 

46223 

TW 

27.2 

473 

4-28-69 

46266 

LC 

27.3 

458 

4-29-69 

46267 

TW 

27.0 

458 

5-  2-69 

46268 

TW 

27.2 

447 

5-  2-69 

46268 

FA 

27,8 

469 

5-  2-69 

46269 

TW 

27.0 

503 

5-  6-69 

46^  0 

TW 

26.7 

490 

5-  9-69 

4627*1 

LC 

27.2 

455 

5-12-69 

46272 

TW 

26.9 

451 

5-12-69 

46273 

LC 

27.6 

450 

5-14-69 

4627*4 

TW 

26.9 

465 

5-16-69 

46275 

LC 

27.3 

460 

5-19-69 

46276 

TW 

27.0 

470 

5-21-69 

46278 

TW 

27.3 

461 

5-22-69 

46279 

TW 

27.5 

455 

5-27*r69 

46280 

TW 

27.7 

439 

5-27-69 

46281 

LC 

28.1 

464 

6-  3-69 

46282 

LC 

27.5 

457 

6-  4-69 

46283 

TW 

27.2 

466 

6-  5-69 

46284 

TW 

27.6 

467 

6-  9-69 

~46285 

LC 

28.0 

453 

p'a$&/  ,  •  _•  .Weight 

Fi&eS'  Lot  No.  Ccenponenfcs  {(Sr&lzis) 


ChWvber  Pres* 
x  180  ^si 


6-10-69 

46286 

TW 

5,7.9 

450 

6-13-69 

46287 

LC 

27.9 

445 

6—16—69 

46288 

LC 

27.4 

455 

6-17-69 

46289 

TW 

27.3 

462 

6-17-69 

46290 

LC 

27.6 

449 

6-23-69 

46337 

LC 

27.7 

448 

6-18-69 

46338 

TW 

27.2 

445 

6-20-69 

46339 

LC 

27.3 

471 

6-24-69 

46340 

TW 

27.4 

481 

6-30-69 

46341 

LC 

27.0 

458 

7-  1-69 

46342 

LC 

27.6 

451 

7-  1-69 

46343 

LC 

27.8 

453 

7-  2-69 

46344 

TW 

27.0 

477 

7—  3—69 

46346 

TW 

27.6 

460 

7-  7-69 

46346 

TW 

27.5 

442 

7-  7-69 

46347 

LC 

27.8 

473 

7-  8-79 

46348 

TW 

27.4 

465 

7-10-69 

46349 

LC 

27.6 

473 

7-14-69 

46350 

LC 

26.9 

466 

7-16-69 

46351 

TW 

27.5 

469 

7-17-69 

46352 

TW 

27.8 

481 

7-23-69 

46353 

LC 

27.9 

498 

7-24-69 

46354 

TW 

28.0 

479 

7-28-69 

46355 

LC 

28.3 

467 

7-28-69 

46356 

TW 

28.2 

469 

8-  5-69 

46357 

TW 

27.3 

482 

8-  8-69 

46358 

LC 

27.8 

467 

8-11-69 

46359 

TW 

27.5 

456 

8-12-69 

46360 

LC 

27.5 

470 

8-13-69 

46361 

TW 

27.6 

462 

?;»*>»■<-  n7a^TV:s’<>W’K*i  g-;<,/IW&w»«r»  *  **»yct»j 


TX3XV7&1  V  v1>  V^UA-'OR^W^  V3<^^ 
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Date^’"^; .  . '  ’■ 

Fired  -Lot’  %>  .* 

.  *  > . 

•  ?  1  >  ’ 

Components 

-  Charge 
Weight 
(Grains)* 

Corrected  Mean 
Chamber  Pressure 
k  100  psi  .  ' 

8-14-69 

46362 

LC 

27.9 

451 

8-15-69 

463*63 

TW 

27.8 

459 

8-18-69 

4  6364’ 

LC 

27.9 

467 

8-19-69 

4636!> 

TW 

27.9 

461 

8-19-69 

463<66' 

LC 

27.9 

468 

8-20-69 

4641;1 

TW 

27.7 

458 

8-21-69 

46412* 

LC 

27.9 

475 

9-  2-69 

46413 

TW 

27.6 

475 

9-  3-69 

46414 

LC 

27.8 

461 

9-  8-69 

46415 

TW 

27.6 

493 

9-10-69 

46416 

LC 

28.0 

466 

9-11-69 

46417 

TW 

28.4 

491 

9-15-69 

464i8 

LC 

27.5 

442 

9-16-69 

46419 

LC 

28.1 

456 

9-19-69 

46420 

LC 

28.6 

478 

9-22-69 

46421 

LC 

28.3 

426 

9-22-69 

46425 

LC 

27.6 

445 

9-23-69 

46424 

TW 

27.8 

455 

9-25-69 

46426 

LC 

27.8 

458 

9-26-69 

46422 

LC 

28.4 

455 

9-26-69 

46423 

LC 

28.3 

453 

9-30-69 

46427 

LC 

27.9 

452 

10-  2-69 

46428 

TW 

27.5 

486 

10-  7-69 

46429 

TW 

77 . 6 

472 

10-10-69 

46430 

LC 

28.3 

477 

10-2G-69 

46422 

LC 

27.9 

462 

10-20-69 

46423 

LC 

28.0 

465 

10-21-69 

46884 

LC 

27.9 

476 

10-22-69 

46885 

LC 

27.8 

465 

10-28-69 

46881 

LC 

28.1 

480 

10-29-69 

46882 

TW 

27.7 

460 

200 


Date  v. 
Pir^d 

•tot  $G. 

.  Cpaponent^s 

Charge 
’  Weight 

V(  (Grains) 

Corrects  Mean 
Chamber  Pressure 
x  100  psi 

t 

r**"*  ■  •  *  $  * 

•:  •  s  r 

i;  *  * 

ti.  —  rt.-it 

10-30-69 

46683 

LC  -  *- 

*  27 19 

458 

11-  7-69 

46806 

LC 

27.6 

459 

11-17-69 

46888 

TW 

27.0 

448 

11-21-69 

46890 

LC 

27.1  . 

447 

12-  8-69 

46891 

LC 

26.9 

452 

12-11-69 

46892 

TW 

27.0 

451 

12-16-69 

46893 

LC 

.  26,8 

452 

12-17-69 

46894 

TW 

26,2 

443 

12-22-69 

46895 

TW 

27.2 

442 

12-29-69 

46896 

TW 

26,6 

M2 

11-18-69 

46889 

TW 

27.3 

463 

1-  5-70 

46897 

TW 

26.5 

mo 

1-  7-70 

46898 

TW 

26.3 

458 

1-  8-70 

46899 

TW 

27.1 

460 

1-13-70 

46900 

LC 

26.9 

457 

1-15-70 

46506 

TW 

27,1 

435 

1-29-70 

46507 

LC 

27.5 

469 

2-11-70 

46508 

LC 

27.3 

453 

2-19-70 

46509 

LC 

26.5 

481 

2-19-70 

46510 

LC 

26.3 

481 

2-20-70 

46511 

LC 

26.5 

461 

2-20-70 

46512 

LC 

26.6 

465 

2-20-70 

46513 

LC 

26.8 

465 

2-24-70 

46514 

LC 

27.6 

464 

2-25-70 

46515 

LC 

27.2 

456 

2-27-70 

46516 

LC 

27.0 

488 

3-  9-70 

46517 

LC 

27.2 

479 

3-11-70 

46518 

LC 

27.1 

462 

3-17-70 

46519 

LC 

27.0 

475 

3-23-70 

46520 

LC 

27.2 

470 

3-26-70 

46521 

TW 

26.7 

440 

nseM  6&io-3'i2r>^  -jp-.s  O 

ssi/'a3si?  .f.30ir»i.*,'w-  jr.g.o't  charge  Corrected  ] 

Date  or*1  7  °’w &mc  'ChknbhrWm 

Fired  ~Eot  So .  Components  “(Grains)  x  100  pel 


4-  1-70 

40522 

TW  i-  .  " 

27,4 

461 

4-13-70 

46523 

LC  . 

27.1 

463 

4-22-70 

46324 

TW  ■  .  •  f. 

27.1 

464 

5-  4-70 

46526 

FC 

27.2 

436 

5-  4—70 

40527 

TW 

27.0 

452 

5-  7-70 

46528 

TW 

27.1 

485 

5-12-70 

46529 

LC  < 

27.8 

-454 

5-22-70 

46530 

TW 

27.3 

448 

6-  2-70 

4&604 

RA 

28.0 

489 

6-  3-70 

46605 

WW 

27.9 

422 

6-  5-70 

45606 

ww 

28.2 

434 

6-17-70 

46607 

LC 

28.0 

.467 

6-18-70 

46608 

WW 

28.3 

449 

6-25-70 

46609 

LC 

27.9 

458 

4-23-70 

46525 

LC 

27.1 

471 

6-26-70 

46610 

TW 

27.7 

455 

7-  7-70 

46612 

TW 

27.4 

452 

7-  8-70 

46611 

LC 

27.4 

470 

7-17-70 

46613 

TW 

27.2 

466 

7-20-70 

46614 

LC 

27.4 

468 

7-20-70 

46615 

LC 

27.3 

466 

7-21-70 

46616 

TW 

27.1 

459 

7-28-70 

46617 

TW 

26.9 

464 

7-29-70 

46618 

TW 

27.1 

457 

7-3C-70 

46619 

LC 

27.5 

452 

o 

r" 

i 

i 

CO 

46620 

TW 

27.3 

465 

8-  4-70 

46621 

TW 

27.8 

460 

8-  7-70 

466 22 

TW 

27.6 

466 

8-13-70 

45623 

&c 

27.8 

492 

8-20-70 

4(5624 

TW 

27.6 

456 

202 


-  '  "*  Charge  Corrected  ftean 

Date  >  v.  Weight  Chamber  Pressure 


Fired 

hot  No. 

Components 

(Grains) 

‘  x  100  psi 

8-24-7 C 

46625 

TW 

27.5 

455 

8-31-70 

46626 

LC 

28.1 

443 

8-31-70 

46627 

WW 

28.1 

437 

9-  2-70 

46628 

LC 

27.4 

464 

9-  2-70 

46938 

TW 

27.4 

468 

9-10-70 

46939 

LC 

27.9 

481 

9-10-70 

46940 

TW 

27.5 

443 

9-14-70 

46941 

LC 

27.7 

500 

9-17-70 

46942 

TW 

27.3 

447 

9-21-70 

46943 

LC 

27.2 

456 

9-23-70 

46944 

TW 

27.6 

465 

9-30-70 

46945 

TW 

27.9 

450 

9-30-70 

46946 

TW 

27.6 

474 

10-  5-70 

46948 

TW 

27.4 

452 

10-  9-70 

46949 

TW 

27.0 

439 

10-  9-70 

46947 

LC 

27.9 

428 

10-12-70 

46950 

RA 

27.4 

445 

10-14-70 

46951 

LC 

27.3 

469 

10-15-70 

46952 

LC 

27.4 

470 

10-22-70 

46953 

TW 

27.1 

445 

10-27-70 

46954 

LC 

27.5 

464 

10-27-70 

46S55 

LC 

27.3 

485 

10-27-70 

46956 

TW 

27.4 

461 

10-29-70 

4695/ 

LC 

27.3 

473 

10-30-70 

4695-3 

TW 

27.2 

463 

11-  5-70 

46559 

WW 

27.2 

4, <4 

11-  9-70 

46960 

TW 

27.1 

452 

11-10-70 

46961 

LC 

27.0 

490 

11-10-70 

46962 

TW 

27.4 

462 

11-12-70 

46999 

LC 

27.0 

481 

203 


SOS 


RiiOK  tyxiC-O1}  ‘It.  ■  <xK.Fi" 

Date/eesaq  ns-draj/f ;  Jai 

Fired  i  i»ptC'  No .  Components 

SKSrge  "Corrected  Mean 
Weight  Chamber  Pressure 
(Grains)  .  •  *  100  .psi  ,  : 

11-17-70 

47000 

TW  .  , 

27.0 

•47-6 

11-17-70 

47,001 

LC 

27.0 

451  i 

11-24-70 

47002 

TW 

27.1 

4Stf 

12-  1-70 

47003 

LC  *  . 

27.0 

-45.5' 

12-  4-70 

47004 

TW  ■ 

26.8  ' 

458' 

12-  8-70 

47005 

LC 

26.0 

489 

12-10-70 

47006 

TW 

26.5 

480 

12-10-70 

47007 

LC 

26.2 

480 

12-15-70 

47008 

TW 

26.5 

466 

12-18-70 

47.009 

TW 

26.3 

464 

12-22-70 

47010 

LC 

26.3 

441 

12-23-70 

47011 

TW 

26.3 

447 

12-23-70 

47012 

LC 

26.5 

467 
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APPENDIX  V 
COMPUTER  PROGRAMS 

In  this  appendix,  different  programs  used  for  calculations  In  this 
report  are  included.  The  programs  are  written  for  UNIVAC  1108  computer. 

TREND 

Purpose:  To  estimate  the  underlying  process  from  a  given  series  of 
observations  using  Semi  Average,  Cumulative  Sum,  and 
Moving  Average  techniques . 

Input: 

K(I)  =  different  values  of  constant 

AK(II)  =  different  values  of  constant  K2 

KX(I)  =  different  values  of  constant  K. 

X  (I)  =  observed  values  of  input  variable 

(e .  g . ,  chamber  pre  s  sure) 

Output:  Plots  of  estimated  process  behavior  based  upon  the  three 


methods  of  estimation. 


TREND 


TREND  DETERMINATION  IN  TIME  SERIES 

DIMENSION  X(1000) .XBAR(IOOO) .Z(IOOO) .SUM(lOOO) .XXBAR(lOOO)  . 
1K^)*AK(5)*<K(5) 

READ  i.N.O  (J).  J  =  1*3) 

1  FORMAT(4I10 ) 

READ  2. ( AK ( I )  I  =  1.3) 

2  FORMAT ( 3F10«0 ) 

READ  3.  ( KK < I  )  *  1=  1.3) 

3  FORMAT (31 10) 

READ. ( X( I ) »  I*  l.N) 

DO  4  I  «  l.N 

4  X ( I >  =  X(I)*100. 

Do--#  I  -*JOlVN’ 

5  Z(I)  =  I 
CAkU;^d'NITPL(9dOV8)  - 
METHOD  OF  SEMI  AVERAGES 
DO  100  J  «  1.3 

Y  =  K(J) 

JJ  =  1 
L  =  N/K(J) 

DO  6  l.L 

SX  =  0. 

SS  =  0. 

JJJ  =  JJ+K  { J )  *—  1 
DO  7  I  =  JJ.JJJ 

7  SX  =  SX+X(I) 

X  BAR ( M )  =  SX/Y 

6  JJ  =  J  J+K  (  j  ) 

CALL  GRAPH ( Z  »3HL IN»'XBAR.3HLIN.L  .4HN0NE .5HS0L ID. 

1 2  HS>$  »2H$$  .4HAXES  »2H$$ .4HFULL .4HNULL ) 

100  CONTINJE 

CUSUM  CHART 
DO  200  J  =  1.3 
S  =  0. 

DO  8  I  s  l.N 
S  =  S+X(I)-AK(J) 

8  SUM ( I )  =  S 

CALL  GRAPH ( Z .3HLIN .SUM.3HLIN  »N .4HN0NE . 5HS0L ID. 

12H$S,2H$$ .4HAXES  »2H$$ .4HFULL .4HNULL ) 

200  CONTINUE 

METHOD  OF  MOVING  AVERAGES 
DO  300  J  =  1.3 
YY  =  KK(J) 

KKK  =  KMJ) 

SS  =  0* 

DO  9  I  =  l.KKK 

9  SS  =  SS+X ( I ) 

XX BAR  (  .1 )  =  SS/YY 
NN  =  N-KK(J)+1 
DO  10  I  =  2 .NN 
LL  =  I +KK ( J) — 1 
SS  =  0. 

DO  11  M  «  I .LL 
11  SS  =  SS+X ( M ) 

XX BAR ( I )  a  SS/YY 

10  CONTINUE 

CALL  GRAPH ( Z  » 3HL IN.XXBAR  »  3HL IN  »NN .4HN0NE  »  5HS0LI D » 

12H$$  »2H$$,4HAXES  »2H$$ .4HFULL .4HNULL ) 

300  CONTINUE 

CALL  ENDPLT 

END  _ 
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HISTOGRAM  AND  CUMULATIVE  DISTRIBUTION  FUNCTION 

Purpose:  To  obtain  one  dimension  histogram  and  cumulative  distribution 
function. 

Input: 

Z  (I)  =  observed  values  of  input  variable 

(e.g.  Chamber  Pressure) 

N  =  Number  of  values  of  input  variable 

JX  =  Number  of  groups  in  which  the  variable  range  is 

divided  to  obtain  the  histogram . 

Output:  Plots  of  histogram  and  cumulative  distribution  function 
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•  F OR  *  S I  HliSTO 

SUBROUTINE. HI$TO(N*Z) 

DIMENSION! (5000).Y(50).X(50)»AY(50)»BY(50) 

DIMENSION  FEU  16)  **  FELCH 

DATA  FEL  /  6HHIST0G»6HRAM  OF.6H  0BSER*6HVATI0N.6HS$S  ♦  **  FELCH 

16H0BStRV,6HATl0NS,6H$S  .6HNOMBER.6H  OF  OC .6HCU&ENC »6HESSS  *  **  FELCH 

26HHIST0G,6I  *AM  OF.6H  RESID *6HUALS$S  /  **  FELCH 

N6=i  -  }'■  .  .  . 

JX  =  20 

CALL  URSRCH<O.N»Z.IBIG,B.O»DUMMY)  **  FELCH 

CALL  URSRC.H]  1  *Z *^SMA?iS *0  .DUMMY  1  **  FELCH 


RANG=Z( IBIG)-Z(ISMA) 

Q**  4#*****#**#*##*####*#tf###*#***#**##*****###**######****##*tf 

PRINT  50»Z(IBIG)»Z<ISMA) 

50  FORMAT ( 10F10.3 ) 

AJX=JX 

del=rang/ajx 

X ( 1 ) =Z ( ISMA ) 

C*#  5## +####*#■*###*##*•*#*■»*•#######■»«•  *■#**■**•»•<)•**■&##*#*#■*#**####** 

YU)--  O.OOOUOOOOE+OO 
JX  =  Jx+l 

IF(JX.LT.l)  GO  TO  90 
DO  80  I  =  2  * JX 
Y ( I ) =  O.OOOUOOOOE+OO 
70  X(I)=X(  I-D+DEL 
80  CONTINUE 

X(JX)  =  Z ( 1 3 1 G ) 

90  CONTINUE 

Q**  #****#**#*#  #***#**#*#**#  *######  ***##*#*■##**!■#*&«###*■«•*»*■& 

PRINT  50.  RANG. DEL .X ( 1 )  «X(JX) 

I F (  N  .LT.  N6)  GO  TO  160 
DO  150  I*  N6.N 
I F (  JX  .LT.l)  GO  TO  130 
DO  123  K  =  2 ♦ JX 
I F ( X ( K ) )  lio,  110.  100 
100  IF (Z ( I )-X ( K) )  140.  140,  120 
110  IF (Z ( I )-X ( K) )  140.  140,  120 
120  CONTINUE 
130  CONTINUE 
GO  TO  150 

140  Y(K)=Y(K)+  l.OOCOOOOOE+OO 
150  CONTINUE 
160  CONTINUE 

C*****#*H0LLERITH*C0NSTANT*L0NGER*THAN*6*CHARACTERS************* 

PP  =  0, 

IF ( PP ) 27 » 2*  .28 

27  CONTINUE 

CALL  INITPLt  9,10.8 ) 

C*200  CALL  GRAPH 1 ( X ♦ Y ♦ JX ♦ 6H0X6  .6HAUT0  .32HHISTOGRAM  OF  OBSERVATIONS. 

C*  1.  ♦ 16H0BSERVATI ONS. •  .24HNUMBER  OF  OCCURENCES..  ) 

200  CALL  GRAPHIX.l , Y » 1 » JX .4HN0NE , 5HS0LID ,FEL C 6 > ,FEL(9) .1 »FEL( 1 ) ,7.5rl0*#  hELCH 
l.U)  **  FELCH 

28  CONTINUE 
230  FORMAT (1H1) 

290  FORMAT (10X»F7«2» 1X.3HT0  »F7.2 .1 IX »F6.0»/ ) 


BY ( 1 )  =  AY ( 1  ) 

DO  2021  =  2  *  JX 

202  BY(I')  s'  BY  (  I-li+AY  (T)’:  ’•  ' 

CALL  GRAPH (  X  ,3HL IN  ,BY , 3HL IN » JX ,4HN0NE ♦ 5HS0LI D ♦ 
1'0BS£^VATJ.0NS.SS»  ♦•CUMULATIVE  PROBABI L I TY$$,‘ ♦ 

14HAXES,  ’EMPERICAL  CUMULATIVE  DISTRIBUTION  FUNCT  ION  j>i> '  ,7*5, 
11U.0) 

CALL  ENDPLT 
I F (  JX  *LT* 1 )  GO  TO  340 
DO  330  1  =  1, JX 

C*#  10*#*#*** **##«*■***  ###■##  -tHHHMUHf  *####**■  #*-tHt-*#*#**#  #**##***##** 

310  PRINT  290  *  X  ( I  >  ,X(  1  +  1)  ,Y(  I  ) 

32C  CONTINUE 
330  CONTINUE 

C#*  1  1  *#******##■**#*■&*■«•■«■  ##****#****#•»■##•»*#*#***■«■#*■>(■■***#  #****#•«•#■»• 

PRINT  400 

400  FORMAT!  10X,*1  IMPROBABILITY  ,10X,22HCUMULAT  I  VE  PROBABILITY,/  /  ) 

DO  402  I  =  l.JX 

PRINT  401 ,AY ( I ) »  BY ( I ) 

401  FORMAT (8X,Fl0o6»16X,F10«6) 

402  CONTINUE 
340  PRINT  230 

RETURN 

END 
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BIVARIATE  HISTOGRAM 

Purpose:  The  program  generates  a  two  dimensional  frequency  table  of 

.  i  ,  1 

chamber  pressure  and  port  pressure  values.  This  table  is  a 

■  *  *  • 

representation  of  the  two  dimensional  histogram. 

.  i 

.  •  1 

Input:  (  > 

N  =  Number  of  chamber  pressure  values  also  equal  to  the 
number  of  port  pressure  values  i 

K  =>  Number  of  groups  in  which  the  chamber  pressure  range 

5 

is  divided 

M  =  Number  of  groups  in  which  the  port  pressure  range 
is  divided 

» 

XMIN  =  Minimum  Chamber  pressure  value 
XMAX  =  Maximum  chamber  pressure  value  , 

i 

XXMIN  =  Minimum  port  pressure  value ! 

XXMAX  =>  Maximum  port  pressure  value 

i 

D  “  One  standard  deviation  of  chamber  pressure 

!  ,  I 

DD  =  Half  a  standard  deviation  of  chamber  pressure 
Z(I)  -  Chamber  pressure  values 
ZZ(I)  =  Port  pressure  values 

Out*  Jt: 


Two  dimensional  table,  where  the  entire  set  of  chamber  pressure 
and  port  pressure  values  are  grouped  into  groups  one  standard 


i 


r 

f  j  ♦ ,  i .  -  ,  ' '  * . v  ■  .  -  . 

i  . 

'  I  ’  '  4  • 

deviation  wide  aldng  the  chamber  pressure  axis  and 

»  ,  .  ’ 

half  a  standard  deviation  wide  along  the  port  pressure 

•  i  i 

axis. 

{  i  . 
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C  Qf5|IVARlAT£  HISTOGRAM 

DIMENSION  Z(1000)>2Z(1000)  »Y( 50*50 )  »X  150)  »XX(50J 
READ  20»N*K *M 
20  FORMAT (3110) 

READ  30*XMIN#XMAX*XXMIN*XXMAX»P*DD 
30  FORMAT ( 6F10®0 ) 

READ  1»(Z(I)*  I«1*N) 

1  FORMAT ( 7F10*0 ).,  v  ,  re  ■„  -  .. 

READ  2»  (22 (I.) »  I  *  1*N> 

2  FORMAT (7F10..0.J, 

DO  9  1  =  1  j  N 

9  22(1)  =  .1^22 ( I ) 

DO  5M  »  1«M  ■  ** 

DO  51  J  =»  1*K 
51  Y(I.J)  =  0# 

50  CONTINUE 

x ( i )  »  xMir 
X ( K )  *  XMAr 
X X ( 1 )  a  XXMIN 
XX ( M )  =  XXMAX 
KK  =  K-l 
DO  3  I  =  2»KK 
3  X ( I )  =  X ( 1-1 ) +D 
MM  =  M-l 
DO  15  I  *  2tMM 
15  XX ( I )  =  XX ( I — 1 ) +DD 
DO  4  I  =  2  »M 
DO  5  II  =  ltN 

IF(ZZ( 1 1 )#LE.XX( I )«AND.22( 1 1 )*GT*XX( 1-1 ) )  60  TO  6 

GO  TO  5 

6  DO  7  J  *  2,K 
IF(2( II )-X{J) )  8  »8  ♦  7 

7  CONTINUE 

8  Y ( I » J)  =  Y(  ItJ)+l. 

5  CONTINUE 

4  CONTINUE 

DO  10  J  =  2»K 

PRINT  1 1 » ( Y ( I »J) *  I  ~  2  »M ) 

11  FORMAT (2X»24F4*0) 

10  CONTINUE 
END 


b 

IK* 


0- 
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AUTO  AND  PARTIAL  CORRELATIONS  (IDENTIFICATION) 


Purpose:  For  a  given  series  of  data,  the  program  calculates  auto  correlation 
and  partial  correlation  functions.  These  functions  are  used  to 
identify  the  form  of  the  tentative  time  series  model.  The  same 
program  is  also  used  in  diagnostic  checking  of  residuals . 

Input: 

N  =  Number  of  data  values 
Z(l)  =  One  dimension  array  of  ordered  observations 
KK  =  Number  of  auto  correlations  or  partial  correlations 
required  +  1 

Output*  Plots  of  auto  correlation  function  and  partial  correlation  function 
with  the  respective  significance  (2  sigma)  limits. 


102 

103 

C 


SUBROUTINE  AUTO<N*7.»KK»AR01»AR02) 

,-D I  MENS  I  ON'Z  (1000  >  »RPLOT( 1COO  )  *  IPLOT  (1000  )  »C  (  }  >)  »R  ( 99  >  *  SCUMS  (  2  )  - 
ISCPA'RSt  10)  »T(50*50)»VAR(9y )  *STSR  (  99  ) » STSR 1  ( 99  i  *STSR2 ( 99 ) ♦ I PLQTV ( 
♦TPLOT  (  50  T»VART  ( 5? )  »i>TSP(  5  J )  *STSP1  (  5  J  )  *STSP2(  50 )  *U(  50 )  * S  ( 99  ) 

KKFlX=Kk 


X-K=KK-1 

2BAR=0. 

XN=N-  •  .  . 

DO  '  102  I =1 »N  • 

Z3AR=Z3AR+2( I ) 
ZDAR=Z3AR/XNi 
C0=0. 

DO  103  .1  =1  »N 

co=c6+ <ZTT ) -zsar ) ** 2 


CO=CO/XN 

calculation  OR  R-- 

RPLOT ( 1 ) =0» 


RPLOT ( 2 ) =1 •  , 


RPLOT ( 3 ) =  0 • 

IPLOT ( 1 ) =0 
IPLOT (21=0 
I  PLOT  <  3 ) =C 
DO  104  K-l »KK 
C  C  K )  =  0 . 

NN=N-K 

DO  105  J=1 *NN 

.105  C(K)=C(K)  +  (Z(  J)-ZBAR)*(Z(  J+K>-ZHAR) 

C(K)=C(K)/XN  . 

R ( K ) =C ( K ) /CO 
J=3*K+1 
IPLOT (  J )  =K 
RPLOT ( J ) =0 • 

J=J+1 

IPLOT (J)=K 
RPLOT ( J ) =R ( K ) 

J  =  J  +  1 

I  PLOT { J ) =K 
RPLOT ( J) =0« 

104  CONTINUE 
AR01=R( 1) 

AR02  =  R ( 2 ) 

:<;<LAG=3*XK 
5CLi;-:om=-i. 
oCLIMSI 2 )  =  !• 

C  CALCULATION  of  t 

I F  (  K.K-5C  )1 06*1  06*107 

107  KkX=:o~i 

CO  TO  109 
106  XKX-k.<.-1 

109  T(1.1)=R(1) 

T(2.2)=(R(2)-RC 1 )*»:)/( l.-R(l)*»2) 
T(2»l)=T(l»l)-T(2*2)*T(lfl) 

DO  20  3  ;<=2,XK.< 


A-  0  • 


DO  202  J~1 »K 
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214 


\ 


A=A+T  {  J)  *RI  ;<+l-J)  . 

202  ^=3+T(K»J)*R{  J) 

A=R ( R+l ) -A 

0=1. -D 

T  ( ts  +  1  »K.+1  )  =A/ti 
DO  203  J  =  1  *  K 

203  T(K+1»J)  =  T(K»J)-T(K+1«K+1)#T(K* K-J+l  ) 

C  CALCULATION  OF  VAR  AND  VART 

I PLOTV ( 1 ) =0 
VAR ( 1 ) =1 • /XN 
STSR( 1 )=SORT( VAR( 1 ) ) 

5(1)=R(1)/3T5R(I) 

5TSR1 ( 1 ) =1 • ?6*STSR { 1 ) 

5T5R2 ( 1 :=-l.?6*STSR( 1 ) 

A  =  2  o /XN 
DO  204  .<=2,(U 

VAR(K.)=VAR(K-i>+Av(R(K-l)**2  > 

IPLOTVK>=K-l 
3T  SR ( K )  =3QRT  (  VAR  ( '< )  ) 

STSRl(lO  =  1.96»ST*R(K) 

3T3R2 ( K  5  =-S  T3R1 ( K ) 

S(K)=R(K)/3T3R(K> 

204  CONTINUE 
kk;<=k;«+i 
tplotj 1 ) =0. 

TPLOT ( 2 )  =  1  • 

TPLOT { 2 ) =0. 

DO  20 '3  Kal.SU/. 

A*1 •/(«-<) 

VART(K)=A 
OTSPt K) =3GRT « A) 

STGPl(K)al.96»ST^P(K) 

3T3P2 ( K) =-ST 3P1 (K  ) 

J  (  K  )  =  T  ( ;< ,  K  )  / 3  T  3P  (  K  ! 

J  =  3  *  K  +  l 
TPLOT ( J ) =  J . 

J-J  +  l 

TPLOTt J)=T(K.<) 

J  =  J  +  1 

TPLOT (Ji=0. 

203  CONTINUE 

RLAo=  3 '< .< K 
C  PRINT  GUT 
PRINT  300 

f-  OR  MAT  (  1M1///1v.a  ,3SH  cA.lPLb  AUTOCORRELATION  COcF c IC I tNTS// ) 
PRINT  JO  I 

3  v  I  FORMAT  l  2  OX  *  1  7H  AUTO-COcFf  I  C  !  EN  T  *5Xr2^H  Of*  I  F  !  r  0  COEFFICIENT 

PRINT  302*  { I »R(  I  )  »f ( I ) ,  1=1 »\K  ) 

302  r  OR  MAT  i  1  3X » 1  J  » F  1 5 . 4  ^x  »  F 1 3  «  4  ) 

PRINT  303,  ZUAN.CO 

303  FORMAT (  10X»19h  MEAN  OF  THE  aERlti  »F14.4/  IOX, 

123H  VARIANCE  CF  THE  oLRIES  .F1C.4) 

PRINT  304 

0*4  FORMAT  (1h:///1vm,27H  OA^PLc  PARTIAL  COR  RtL  A  T  I  ON  //) 

PRINT  302 

j'J  'j  r  J  ;  T 1 A  T  (  2  jX  *  1  7  HP  t\ k  T  -C'JttRi  L  M  1  o  .  »  3  X  *  2j  n  u  N 1  f  1 1 J  COcFrlCicNT/) 


/) 


A  fC 

PRINT  302 »  (It  T(Itl)*  U(I)t  I*ltKK<) 

KK=KKFIX 

RETURN 
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US 


ESTIMATION  AND  DIAGNOSTIC  CHECKING 

.  *' '*  -  ~  ~  ,1  nijo:,inon'  P  '•)*•*  1  -  ‘  ( 

von;  H  io  .  u.sv'oq  jrfi  m  xio*‘jq2 

'll  ,  '  O *>'**!  i  c  V 

Purpose:  The  prog ram  estimates  the  parameters  in  the  proposed 

•  time  series  model  by  the  method  of  nonlinear  least 

io  .  ■  -jfu  ■<  u,  •  *- ■  ...  -  ,o 

squares  ;  ,  Thertestofor  the  adequacy  of  the  model  is 
done  by  considering  the  autocorrelation  and  the  par- 

j  N 

tial  correlation  functions  of  the  series  of  residuals, 

’  i  ' "j 0  •  .  .  ,vl  -  ;•  ",  - 

and  a  chi-square  statistic  based  on  the  autocorrela¬ 
tion  function. 

f  t  *  -  , 

Input : 

A  general  time  series  model  can  be  written  as 
(l-<t>1B-<..<frpBP)  (1-^^'B-;.  )  (l-BS)dl(l-B)d(Zt-y) 

-  =6o+(l-01B-...-0qBq) (lr01,B-...-0qiBql)at 

NREP  =  Number  of  models  to  be  fitted 

NDR  =  Number  of  observations  in  the  data  series 

Z (I)  =  One  dimensional  array  of  data  values 

MAXI  =  p  +  p'  +sdx  +  d 

NP  =  Number  of  parameters  to  be  e<  timated 

NRD  »  d 

NSD  =  s 

NSEA  = 

INC (J)  =  One  dimensional  array  of  size  6.  It 

contains  the  information  regarding  the 
number  cf  parameters  of  each  type  as 
specified  by  the  model. 


IOPA  ( J)*~«  Ons  dimensional  array  of  size  NP.  It 
specif ieo  the  powers  of  B  associated 
with  each  parameter  in  the  model,  in 
3u •vqo-'cY  *  ;  vi  r  the  sequence  from- left  to  right. 

'  .iPA(J>  a  one ^dimensional  array  of.  size  NP, 

specifying  the  initial  estimates  of 
s'  '  i*  ■  the  parameters  in. sequence. 


Special  Subroutine  Used: 

s‘  K  -C  •  ^  ’  * 

The  subroutine  UWHAUS  supplied  by  the  University 

,  *  t 

of  Wisconsin  Computing  Center  is  used  for  the 
nonlinear  least  squares  estimation  of  the  parameters. 


Output: 

The  final  estimates  of  the  parameters,  95%  confi¬ 
dence  limits  on  the  parameter  (on  linear  hypothesis) 
and  the  correlation  matrix  of  the  parameters  are 
pointed.  The  program  also  prints  the  autocorrelations, 
the  partial  correlations  of  the  series  of  residuals 
and  a  chi-square  statistic  based  on  the  autocorrela¬ 
tions  , 
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eis 


I 

c 

s 

? 

in 

II 

l* 

?n 


•  U  A  T  I  *  *  * 

#  «  J  t  *  i  f  S,  t!t>« 


v'  » #  O  ^  *  f  '  i  Of  A '«  <  *  V  , 


i  r»fj 


SM' dfeo'ro  k&»  ^  ‘  ».  «‘J. 

,'Cw  .v,fr  -r-  ,. 


v  'biM'EN^I?)k>fAi<9i^)V2’1«)9^)  *PA(*>b)  ,T0PAfS6}  ,!NC(ft)*Nt(th')*YnQO)  *ZP(20 
XO*lO),CLL(?0O,lO>  *£1)1,.  (2C0*  10)*DIFF(50)  ,$IGNS(  SO)  ,RHO<S0',l0)  »ST£l  50 
X,10)  ,F(10)  ♦SCRATC'U'obh)  .WORT  “  ’  '**  '  ‘  *  *  '  * 

COMMON  /TWO/NDR,NRO,NSD,NSFA 

NF=0  J  ' 


NTOsO 

NLOG=0 

nad=o 


ND=NAO 

M0F=10 

mp=io  ,  1  ■’ 

NY=n 

*it=ao 
fpsi=,ooi 
fps?=,.oov  ,  , 

FLAM=‘.01  ‘  -  . 

FNU  =  m. 

NL=m 

RFAD  1*  NRFP,NDR 

RFAD  7,  (7(1)*  1=1 ,NDR)  '  ’ 

FORMAT  (7014) 

DO  70  NN=1*NRFP 

RFAD  ]«  MAXI *NP 

PFAD  1*  NRD*NSD*NSCA 

(1-B)##NRD  *  (  1-B**m<;FA)**M.S0 

IF  (MV  .F.O.  0)  GO  TO  «i 

READ  7*  ( Y ( I ) *  1=1 ,NY  ) 

CON T I  MU F 

RFAD  1 ♦  ( I  NO ( I)  ♦  1  =  1 ,A> 

RFAD  1  *  ( IOPA( I ) ♦  1  =  1 »NP) 

FORMAT ( 7F1 0*0 ) 

RFAD  7*  ( P A ( I ) »  I=1,NP) 

IF  (NTO  .fq,  O)  GO  TO  in 
©FAD  1*  ( NT ( I ) »  1=1, NTO) 

CONTTNUF 

nob=ndr-maxi 

IF  (NL  • FO*  0)  GO  TO  1 A 
DO  11  I  =  1  »NP 
GIGNS( I)=n. 

D I FF ( I ) =0. OOA 
CONTINUE 

NOTMSs(i*MP+?*NP#*?+?*NOR+NOR«NP 

npror=nn 

call  FSTIM(NPR0R* A, NOR, Z,P A, DI^F, SIGNS *FPS1,FPS2»MIT»FLAM*FNU»SCRA 
XTC»NDIMS»NLOG» I NC » IOPA  *NRD,NSD , NSEA ,NL , NAD ,RH0,STE ,E ,MCc  *MP ) 
CONTINUE 

IF  (NF  .FO.  0)  GO  TO  70 
OONTINUF 


FND 


C 


SUBROUTINE  1 
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SU8R0UT INF  PST IM (NPROB *B .NOB.YYiPA *OIFF .SIGNS.EPS1 #EPS2*MIT . FLAM.F 
XNU»SCRATC»NOIMS.NLOG.IC*!P.NA*NB»NC*NL»NAD.RHO.$TF*F»MCF»MP> 
DIMENSION  A(RR9)  »Z(RR9).C<R9G)*Y(OQ<5)  ,PA  { 50 )  •  I0PA{  SO )  ♦  INC  (6 )  »DI  FF( 
X50)  tSlGNSJ  SO)  *RH0(NL»10)  *3TE ( NL»10)  *E  ( 10 )  *$CPATC  (NDIMS  I » B  { 999)  *YY ( 
X090J  *IPj50)»IC,(A) 

COMMON  A,Z.NDR»IOOA,INC*MAX11 .mbo*nrd*nsd.nsea 

pxtfrnal  tsmod 

KK*0 

nrd»na 

NSO=NR 

NSPA«NC 

no  a  j*i *a 

INO( J)=IC( J) 

KK=KK+INC( J) 

IF< I NC ( J )  .LT.  0)  00  TO  AS 

A  continuf 

! F ( NSCA  ,LT.  0)  GO  TO  55 
MM=n 

IF( ( I NC ( ? )  .NF.  0  .OR.  NSD  .NE.  0)  .OR.  INC(6>  .NE.  0)MM=1 
r  TF(MM  .FO.  1  .AND.  NAFA  ,Lr.  1  )G0  TO  AS 

IF(NPO  .LT.  0  .OR.  NSD  .LT.  0)  GO  TO  55 
IF(INC(A)  .GT  «  1  .OR.  INCU)  .GT.  1)  GO  TO  AS 
IF(NRD  .GT.  0  .AND.  INHA)  .FQ.  1)  GO  TO  55 
!F(NSD  .GT.  0  .AND.  INC(A)  .FQ.  1)  GO  TO  55 
NP=KK 

IF  ( NP  .LF.  0  .OR.  NP  .GT.  50)  GO  TO  55 
DO  5  J=1 *nd 
IOPA( J)=IP( J) 

IF(TODA(J)  .LT.  0)  GO  TO  AA 
A  FONT  I  Ml  IF 

IF( NL  .LF.  0  .OR.  NOR  ,LP.  0)  GO  TO  55 
I F { NAD  .LT.  0  oOR.  NAD  .GT.  10)  GO  TO  55 
MAX  1  aNRD'fNSD#NSF  A 
KK=INC(H 

IF (  TNC(l)  .NF.  0)MAX1*MAX1-MODA(IOO 
KK=KK+INC<?) 

IF t INC ( ? )  .NE.  0)MAX1*MAX1+I0PA(KK) 

NORcNOR+MAXl 

IF ( NDR  .GT.  RRR)  GO  TO  55 
XK=KK+INC(A)  +  TNC(A)4-TNC(  A) 

MAX 

I F ( I NC ( 5  )  .NF.  0)WAX?aI0pA(XK) 

IF ( INC( 6 )  .NF.  0)MAX?«MAX2+IOPA{NP) 
mro»maxi 

IF ( MAX?  ,GT .  MAXI )MRO»MAX? 

IFINDR  .LF.  MBO  .OR.  MBO  .GT.  100)  GO  TO  55 
MAXllsMAXl+1 
DO  A  J«1»N0P 
7A  J)aYY(  J) 

A  IFINLOG  .NF.  0)  Z{ J)=LOG(Z( J) ) 

DO  10  J«MAX1I .NDR 
< *J-MAX 1 
m  Y(K)"7(J) 

no  1?  J-1.MAX1 
AiJ)*0 


1? 


220 


TALL  IJWHAUSI NPROB »TSMOD»NOR* Y  *NP» PA*DI FF* S TGNS »EPS1 * EPS2 *MIT *FLAM» 

xfnu.scra  rc ) 
no  14  I=MAX11.NDR 
J= I -MAX  1 

14  CU)aA(I) 

PRINT  17 

17  FORMAT (1H1 .43X.2SHRFSTDUAL  AUTOCORRELATIONS//) 

CALL  ACOR ( C  » NL .NOB .M AD  »RHO *  STF  *  F* MCE  *MP ) 

SlJM=0 

00  ?0  I=MAX1 1 »NOR 
20  SlJM=SUM+(  A(  I)-FC  1 )  >**? 

DT\/=NOR*<  NOR-NP  ) 

VARAR=SUV/OIV 

TNSTO=ARS(FU)  )/SORT(VABAR) 

PPTNT  23  »  TNSTO 

23  F0RMAT(///10X.40HMEAN  OF  ORIGINAL  SERIES  OF  RESIDUALS  IS  »F6«3*31H 
X  STANDARD  DEVIATION*  FPOM  7FRO.) 

DO  J=1 »NOP 
26  P ( J ) = A ( J ) 

GO  TO  60 
SS  OPINT  S3 

S3  FORMAT!  1H1  »10X*24HOA(?AMFTFP  rRROP  IN  FST  I M  )  ■ 

60  RFTUPN 
PNO 

f  S( IBROIIT  INF  2 

SI  1RR0I  IT  I  NF  TSMOD  (  NPPOR  »P A ♦  F  *  NOP  *NP  ) 

DIMENSION  PA(S0)»F(9Q<?),A(999)  » ?.  (  999 ) » I OPA  ( 50 )  *  I NC  ( 6  >  ♦  T 1 100  )  ♦  C  ( 1 00 
X)*CF(100) *D< 10*10) *DS(10*100)*W<999) 

COMMON  A *7  *  NOR ♦ TOP A  *  INC  »MAX1 1 ♦ va0*NRD*NSO»NSFA 
COMMON  /ONF/C,OP 
no  3  .1=1  *MRO 
T  (  J  )  =  0 
C( J)=0 

3  CF(J)=0 
L*o 

IFIMRn  .FO.  0)  GO  rn  )7 
L=NPD 

n( l  ,1 )=-l.o 

I F ( NPn  .FO.  1 )  GO  TO  3 
no  4  j=2.Npn 

4  P( J*1 ) =0  C J-l .1 )-1.0 
3  no  13  j=i,NRn 

no  13  <=2»NPO 
I F ( X  .GT.  ) )  GO  TO  10 

o(  j,K)=r><  j-i  .if)-0(  j-i  .<-i ) 

GO  TO  1 3 
lO  D(J*K)=0 
13  CONTINUE 

00  IS  K=1,NP0 

15  C(X  )=n(NRO.K ) 

17  IF(NSn  .FO.  0)  GO  TO  A  3 
maX*NSO*NSFA 
no  20  j=i»Nsn 
no  70  K  =  1  .MAX 
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?0  0S(j*KJ=0  * 

osn*NSEA)»-i.o 

MlN=?#NSFA 

IF(NSO  .FO.  1)  GO  TO  ?fl 
DO  ?4  j*2*nio 

74  OS ( J* NSEA ) »DS ( J-l  «NSEA)-l*0 
?fl  00  H  J*1 *NSO 

oo  11  KjsMTNsMAX  *NSFA 
I F (  ( K/NSFA 1  .GT.  J)  GO  TO  10 
0S( Jt<)»OS( J-1,K)-0A{ J-l.X-NSEA) 

GO  TO  11 

10  0S(J*K)=0 
1?  rONTTNUF 

l=msfa#nso 
00  17  M=1  ,NRO 
T ( M ) eT ( M ) +C ( M l 

00  17  j*nsea*l*nsea 
I F ( M  .EO.  1 )T<J)*T(J)+DS(NSD*J) 

N= J-4-M 

17  T(N)rT(N)+C(M)*DS(NSO.J) 

L2L+NPD 
OO  40  J=1  »L 
C( J)=T( J) 

40  T  (  J )  =0 

41  MIN=1 
MAX=0 

00  6 0  1=1*6 

TF(TNrU)  .FO.  0)  GO  TO  60 
MAX=MAX+INC<  I ) 

I F C I  .EO.  1  .OR.  I  .FO.  4)  GO  TO  60 

00  48  M=MIN*MAX 

K«IOPA(M) 

T ( K ) =T ( <) —PA ( M ) 

IF(l  .EO.  0)  GO  TO  48 
00  41  J=l*l 

IF ( m  *E0.  MIN  .AND.  T  .LE.  ? ) T C J ) «T C J )+C< J) 
JFiW  .EO.  MIN  .AND.  I  .hQ.  6 ) T ( J) »T ( J  )>CF ( J) 
N»J+IOPA(M) 

IF ( I  .EO.  6>T(N  )*T(NJ-PA(M)#CF|J) 

41  I F <  T  .LE.  ?)T(N)«T(N)«PA(M)#r( J) 

40  CONTINUE 

L»L+IOPA(MAX) 

00  11  J*l*L 

IFU  .GF.  1 ) CF  {  J }  =T  {  J ) 

T F C I  .LF.  ?)C(J)«T{ J) 

11  T(J)»0 

TFi I  .FO.  2 )L»0 
AO  MI^aMIN+TNOI I) 

00  66  J»1«MR0 
CF( J)»-fF( J» 

66  rui».rui 

KK*INrm  +  INCCM  +  lN'*(l> 

00  70  J»1*N0R 

70  IF( INC ( 1 )  .FO.  1 )7( J)a?{ J)-PA<KK ) 

KK°KX+ I NC ( 4 ) 


.  w  V 

IF( !N0(4)  .FG.  DCONSTaPAdOO 

IF(INC(4)  »F0.  (1)C0NST=0 

DO  79  KsMAXII.NOR 

A(K 1=7 (<) -CONST 

!F(<  .GT.  WR'o)  GO  TO  7A 

KK=K-1 

DO  74  J=1  *KK 

74  A  t  K ) =A ( K ) ~C ( J ) #Z ( K-J ) +CF ( J ) #A ( K-J 1 
GO  TO  79 

7 ft  DO  7«  J«l»MRO  •' 

7R  A(K)aA<K>-C(  J>#?(K-J)+CF(  J)*A(k;-j) 

79  CONTINUE 

00  fil  J«1,N09 
81  7  ( J )  =W ( J ) 

MAX  1 =MAX 1 1  —  1 
no  RS  J- 1 « NOR 
LL=J+MAX1 

A?  F  ( J )  =7  ( LL  )  -A  ( LL  ) 

0FTUPN 

FNO 

C  GURROUTINF  1 

SUBROUTINE  ACOR  (Z*KK*NfND*R»STSR*E»MCE*MP) 

DIMFNSION  Z(999)»C(101)tR{l0l) »T( 50*50) *VAR(101) tSTSRI 101 ) ♦ 
1VART(S0)*STSP(S0)  *U(50)*SU01  )  *F(5) 

COMMON  /TWO/NDR*NRD*MSD»NSFA 
7RAR=0. 

XN=N 

DO  TO?  1*1  * M 
10?  ZBAR  =  7RAR+Z(I  ) 

7 RAP=7BAR /XM 
F(1 WRAP 
CO  =  0o 

DO  10?  1=1  *N 

10?  CO=CO+(Z( I )-7PAP)*»? 

C  calculation  OF  R 

no  104  K  =  1 

r<x  )=  o. 

NNaM-V' 

no  in?  j=i*mn 

10S  C(X)=C(K)  +  (7(  J)-ZRAR)-»(Z(  J+tf)-7BAR) 

R(K)oC(K)/co 
104  CONTINUF 

C  calculation  OF  T 

IF( KK-50)104»106*1o7 
107  XKK=S0-1 
GO  TO  109 
104  KKK=XK- 1 
109  T  (  1  *  1  )  =R  (  1  ) 

T(?*?)*(PI?)-R(1 )**?)/( l.-Rll )»*?) 

T ( ?  *  1 ) =T ( 1 .1 >-T (?*?)*T(1 *1 ) 
no  ?n?  *=?»*<* 

0=0. 

A  =  0  , 

no  ?o?  j=i *k 


s& 

>  2.23 

A*A+T(K*J)#R(K+1-J)  ,^,A  -m  ,  .  - 

20?  CX,J)*R( J)  !-  ■  • 

A»R(K+1)-A  (  , 

T(K+1*K+1»=A/R  •  **  ’t  •, 

DO  203  J»1*X 

?03  T(K-M*J)*T(K*J>-T{I<+1,K+1  J#T!K#K-J+1)  ’  1 

C  CALCULATION  OF  VAR  AND  VART  !  "  !  ■.*  ,v  -  ' 

VAR(1)*1./XN  1 

STSRtn»SOPT(VAR(D) 

sm«R(n/sTSRm  ..  ..  , 

A*?#/XM  l 

00  ?04  K=?*KK 

VAR ( K ) =VAR ( K-l ) +A#(R ( X-l ) »*? )  ' 

ST SR ( X ) a SORT { VAR ( K ) )  l‘ 

S!*)=R(K)/STSR(K>  !  ! 

? 04  CONTINUE 

KKK=KKK+l  i  ,  ,  i 

00  ?.0*>  <=1  »>f KK 

A=1*/(N-K)  i  i 

VART ( K )  =A  1 

ST.SPf  <  Ja.SORT  ( A)  1  , 

U(K)=T(K*K)  /ST.SP  (  K  ) 

?0S  CONTINUF 

C  PRINT  OUT  1  ■ 

PRINT  300  k  ,  • 

300  FORMAT  ( 1H1 ///l OX  *36H  SAMPLE  AUTOCORRELATION  COEFFIC IENTS// >  1 

PRINT  301  , 

301  FORMAT  (?0X»17H  AUTO-COEFF I C I FNT  >S;k»20H  UNIFIFD  COEFFICIENT  /) 
PRINT  30?*  (I,R(I)*S(/),I=1*KK) 

30?  FORMAT  (  1  SX  ♦  I  ?  ♦  F 1  S.4  /PX  *Fl  S.  A )  1 

CHIS0=0*  ‘ 

DO  400  1=1*30  i 

400  CHIS0=CHIS0+R( I ) »Rf T ) 

XXN=NDR-NRD-NSD*NSFA 

CHTSO=CHISQ*XXN  ‘ 

PRINT  401 ♦  CHISO  i 

401  FORMAT ( // *5X*  *CHI  SQUARE  STATISTIC  3ASED  ON  30  AUTOtORRELAT TONS  = 

X*Fl?oS.//)  1 

co=ro/xN  !- 

PRINT  303,  7R4R,C0  i 

303  FORMAT!  10X*1#H  MF AN  OF  THE  SERIFS  *F14*4/  10X. 

1 ?3H  VARIANCE /OF  THE  SERIFS  »E10.4>  i 

PRINT  304 

304  FORMAT  !1  HI  lj/  !\ OX  ♦  ?7H  SAMPLE  PARTIAL  CORRELATION  //’> 

DRINT  30S 

30*  FORMAT  (  ?0X » 1 7HPART-C0RREL AT  ION  .  ,SX,?0H  UNIFIFD  COEFFICIENT / ) 
PRINT  302  »  ( I  *  T (  I  » I  )  *  U(I)»  I=1»KKK> 

RETURN 

FND 


